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ABSTRACT 



We introduce a new tool - ^ESOPUS: Accurate Equation of State and OPacity Utility Software - for computing the equation of state 
and the Rosseland mean (RM) opacities of matter in the ideal gas phase. Results are given as a function of one pair of state variables, 
(i.e. temperature T in the range 3.2 < log(r) < 4.5, and parameter R = p/(T/ 10 6 K) 3 in the range -8 < log(ft) < 1), 
and arbitrary chemical mixture. The chemistry is presently solved for about 800 species, consisting of almost 300 atomic and 500 
molecular species. The gas opacities account for many continuum and discrete sources, including atomic opacities, molecular ab- 
sorption bands, and collision-induced absorption. Several tests made on ^ESOPUS have proved that the new opacity tool is accurate 
in the results, flexible in the management of the input prescriptions, and agile in terms of computational time requirement. Purpose 
of this work is to greatly expand the public availability of Rosseland mean opacity data in the low-temperature regime. We set up a 
web-interface (http://stev.oapd.inaf.it/aesopusl which enables the user to compute and shortly retrieve RM opacity tables according 
to his/her specific needs, allowing a full degree of freedom in specifying the chemical composition of the gas. As discussed in the 
paper, useful applications may regard, for instance, RM opacities of gas mixtures with i) scaled-solar abundances of metals, choosing 
among various solar mixture compilations available in the literature; ii) varying CNO abundances, suitable for evolutionary models of 
red and asymptotic giant branch stars and massive stars in the Wolf-Rayet stages; iii) various degrees of enhancement in or-elements, 
and C-N, O-Na, and Mg-Al abundance anti-correlations, necessary to properly describe the properties of stars in early-type galaxies 
and Galactic globular clusters; iv) zero-metal abundances appropriate for studies of gas opacity in primordial conditions. 
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1. Introduction 

In a gas under conditions of local thermodynamic al equilibrium 
(LTE) and in the limit of the diffusion approximation (DA), the 
solution to the radiation transfer equation simplifies and the total 
flux of radiation F as a function of radius r is given by: 



Planckian. However, in the outermost layers of a star the pho- 
ton mean-free path may become so long that the DA conditions 
break down, thus invalidating the use of the RM opacity. In these 
circumstances, a straight arithmetic average of the monochro- 
matic absorption coefficient (Eddington 1922), designated with 
k p , Planck mean (PM) opacity: 
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where T is the gas temperature, p denotes the density, B(r, T) 
is the integral of the Planck function over frequency, and the 
relation 
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first introduced by Rosseland (1924), defines the Rosseland 
mean opacity kr(p, T). Being a harmonic average over fre- 
quency, kr emphasises spectral regions of weak absorption, 
across which the energy flux is most efficiently transported. 

Both LTE and DA conditions are usually met in the stel- 
lar interiors, where collisions dominate the thermodynamic 
state of matter, the photon mean free-path is much shorter 
than the typical scale length of the temperature gradient, and 
the Kirchoff's law applies with the source function being the 



may be more suitable to represent the absorption properties of 
the gas in a simplified version of the radiation transport equation 
(e.g. Helling et al. 2000). 

Both RM and PM opacities are frequency-integrated aver- 
ages, so that they only depend on two independent state vari- 
ables, e.g. temperature T and density p (or pressure P), and the 
chemical composition of the gas. 

In stellar evolution models it is common practise to describe 
the absorption properties of matter with the RM opacity formal- 
ism, adopting pre-computed static tables of kr which should en- 
compass a region of the bi-dimensional space T-p wide enough 
to cover all possible values met across the stellar structure dur- 
ing the evolution, from the atmosphere down to the central core. 
The chemical composition is usually specified by a set of abun- 
dances, e.g.: the total metallicity Z, the hydrogen abundance 
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X, and the partitions {X,/Z} of heavy elements in the mixture, 
which depend on the specific case under consideration. Frequent 
choices are assuming solar partitions {Xj/Z} = {X,/Z} , or deriv- 
ing [Xi/Z] from other constraints such as the enhancement in a- 
elements (expressed by the ratio [a/Fe]), or the over-abundances 
in C and O necessary to describe the hydrogen-free chemical 
profile in He-burning regions. 

In the literature several authors have calculated kr(p, T) for 
different combinations of the state variables and chemical com- 
position. Let us limit here to briefly recall the most relevant 
efforts, i.e. those mainly designed for supplying the scientific 
community with extended and continuously updated RM opac- 
ity databases. 

In the high-temperature regime, i.e. 10 4 K % T ^ 10 9 K, cal- 
culations of RM opacities are mainly provided by two indepen- 
dent teams, namely: the Opacity Project (OP) international col- 
laboration coordinated by Seaton (Seaton 2005, and references 
therein); and the Opacity Project at Livermore (OPAL) being 
carried on by Iglesias, Rogers and collaborators (see Iglesias & 
Rogers 1996, and references therein). Both groups have set up 
a free web-access to their RM opacity calculations, via either a 
repository of static tables and/or source routines, or an interac- 
tive web mask where the user can specify the input parameters 
and run the calculations in real time. 

In the low-temperature regime, 10 3 K 5 T ^ 10 4 K, widely- 
used RM opacity tables are those provided by the research group 
of the Wichita State University (Ferguson et al. 2005 and ref- 
erences therein). A web page hosts an archive of static RM 
opacity tables, for both scaled-solar and a-enhanced mixtures, 
which cover a wide range of metallicities including the Z = 
case. It should be acknowledged the large body of work made 
by Kurucz, who provides, via web or CD-ROMs, all neces- 
sary atomic and molecular data as well as FORTRAN codes to 
calculate kr (see Kurucz 1993abc), in the temperature interval 
10 3 K S5 T S= 10 5 K, for scaled-solar and a-enhanced mixtures. 
More recently, Lederer & Aringer (2009) have calculated and 
made available via the VizieR Service a large catalogue of RM 
opacity tables for C- and N- rich compositions, with the purpose 
to supply RM opacity data suitable for the modelling of asymp- 
totic giant branch (AGB) stars. Helling & Lucas (2009) have 
produced a set of gas-phase Rosseland and Planck mean opac- 
ity tables for various metallicities, C/O and N/O ratios. It is due 
mentioning also the recent paper by Sharp & Burrows (2007), 
who provide an exhaustive and useful review on the thermo- 
chemistry, techniques, and databases needed to calculate atomic 
and molecular opacities at low temperatures. 

Despite the undeniable merit of all these works, the pub- 
lic access to low-temperature RM opacities still needs to be 
widened to account for the miscellany of chemical patterns - 
mostly relating to the photosphere of stars - that modern spec- 
troscopy is bringing to our knowledge with an ever-growing 
richness of details, and also to allow the exploration of possi- 
ble opacity changes driven by any hypothetical chemical com- 
position. The peculiar abundance features in the atmospheres of 
AGB stars (e.g. McSaveney et al. 2007, Smith et al. 2002); the a- 
enhanced abundance pattern of stellar populations belonging to 
globular clusters (e.g. Gratton et al. 2004) and elliptical galaxies 
(e.g. Clemens et al. 2006, 2009); the large carbon overabundance 
and other chemical anomalies of the so-called carbon-enhanced 
metal-poor stars in the Galaxy (e.g. Beers & Christlieb 2005); 
the striking C-N, O-Na and Mg-Al abundance anti-correlations 
exhibited by stars in Galactic globular clusters (e.g. Carretta et 
al. 2005); the chemical composition of the primordial gas after 



the Big Bang nucleosynthesis (e.g. Coc et al. 2004): these are a 
few among the most remarkable examples. 

In this framework, purpose of our work is to greatly ex- 
pand the availability ofRM opacity data in the low-temperature 
regime, by offering the scientific community an accurate and 
flexible computational tool, able to deliver RM opacities tables 
on demand, and with a full freedom in the specification of the 
chemical mixture. 

To this aim, we have developed the ^ESOPUS tool (Accurate 
Equation of State and OPacity Utility Software), which consists 
of two fundamental parts: one computes the equation of state 
(EOS) of matter in the gas phase, and the other evaluates the 
total monochromatic coefficient, k(v), as sum of several opac- 
ity sources, and then computes the Rosseland mean. The EOS is 
solved for m 800 chemical species, including neutral atoms, ions, 
and molecules. The RM opacities take into account several true 
(continuum and discrete) absorption and scattering processes. 
An interactive web-interface (http://stev.oapd.inaf.it/aesopus) al- 
lows the user to run ifiSOPUS according to his/her specific re- 
quirements just by setting the input parameters (T - R grid, ref- 
erence solar composition, total metallicity, abundance of each 
chemical species) on the web mask. 

The paper is organised as follows. Section [2] specifies the 
bi-modular structure of yESOPUS. In Sect. 12. 11 we illustrate the 
basic ingredients necessary to set up and solve the equation of 
state. Numerical aspects are detailed in AppendixlAl Section lZ2l 
indicates the opacity sources included in the evaluation of the to- 
tal monochromatic absorption coefficient. The Rosseland mean 
is presented in Sect. 12. 2. Tl with details on the computing-time re- 
quirements provided in Sect. 12.2.21 Complementary information 
on the frequency integration is given in Appendix [B] The for- 
malism introduced to describe the different ways the RM opac- 
ity tables can be arranged, as a function of the state variables and 
chemical composition, is outlined in Sect. [3] In Sect. [4] we anal- 
yse five relevant cases of RM opacity calculations, characterised 
by different chemical patterns, namely: scaled-solar elemental 
abundances (Sect. |4~TT i, varying CNO abundances (Sect. [43]), a- 
enhanced mixtures (Sect. 14.31 ). mixtures with peculiar C-N-O- 
Na-Al-Mg abundances (Sect. 14.41 ), and metal-free compositions 
(Sect. 14. 51 ). Appendix [C] specifies the general scheme adopted 
to construct non-scaled-solar mixtures. Final remarks and in- 
dications of future developments of this work are expressed in 
Sect.0 



2. The /ESOPUS code 

2.1. Equation of state 

The equation of state quantifies the distribution of available par- 
ticles in the unit volume, in the form of neutral and ionised 
atoms, electrons, and molecules. At low temperatures (T £ 
6 000 K) and sufficiently high densities, molecules can form in 
appreciable concentrations so as to dominate the equation of 
state at the coolest temperatures. To this respect a seminal work 
was carried out by Tsuji (1964, 1973) who set up the theoretical 
foundation of most chemistry routines still in use today. 

In our computations the EOS is solved for atoms and 
molecules in the gas phase, under the assumption of an ideal 
gas in both thermodynamic equilibrium (TE) and instantaneous 
chemical equilibrium (ICE). This implies that the abundances of 
the various atomic and molecular species depend only on the lo- 
cal values of temperature and density, regardless of the specific 
mechanisms of interaction among them. 
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Solving a chemical equilibrium problem requires three gen- 
eral steps. First, one must explicitly define the gas system in 
terms of its physical and thermodynamic nature. For example, 
the classical problem in chemical equilibrium computations is 
to calculate the state of a closed system of specified elemental 
composition at fixed temperature T and pressure P. The nature 
of the physical-chemical model determines the set of govern- 
ing equations to be used in computations. The second step is to 
manipulate this original set of equations into a desirable form, 
to reduce the number of unknowns and/or to fulfil the format re- 
quirements of the adopted computation scheme. The third step is 
to solve the remaining simultaneous equations, usually be means 
of iterative techniques (see, for instance, Tsuji 1963). 

Rather than solving sets of equations, the equilibrium com- 
putation can be formulated as an optimisation problem, such as 
solving the so-called classical problem by minimising the calcu- 
lated free energy of the system (Mihalas, Dappen, & Hummer 
1988). An alternative approach, based on the neural network 
technique, has been recently proposed by Asensio Ramos & 
Socas-Navarro (2005). 

In this study we adopt the Newton-Raphson iteration scheme 
to solve the chemical equilibrium problem of a gas mixture with 
assigned chemical composition, pressure P (or density) and tem- 
perature T. The adopted formalism and solution method are de- 
tailed below. 



2.1.1. Equilibrium relations 

Under the ICE approximation, the gas species obey the equilib- 
rium conditions set by the dissociation-recombination and ioni- 
sation processes. Generally speaking, the chemical interactions 
in the gas between species A and B may involve the simple 
dissociation-recombination process 



A + Bt^AB 



(4) 



in which both forward and reverse reactions proceed at the same 
rate. In the above equation A or B may be an atom, molecule, 
ion or electron. Of course one may postulate more complicated 
chemical interactions such as 

AB + CD ±5 AC + BD 



or 



A + B + C ±5 ABC 



but these can ultimately be reduced to Eq. (0), in the forms of 
simple dissociation-recombination reactions, i.e. 

(AB) + (CD) ±5 (ABCD) i=> (AC) + (BD) 

(A + B) + C fc> ABC . 

From statistical mechanics we know that for any species A and 
B in equilibrium with their compound AB (usually a molecule), 
the number densities ha, «b, and uab are related by the Guldberg- 
Waage law of mass action: 



K AB (T) 



n A n B 

tlAB ' 



(5) 



where Kab(T) is the dissociation constant or equilibrium con- 
stant of species AB, which depends only on temperature. It is 
expressed with 



(27TllkT\ 3/2 gint,A(2int,fl 



K AB (T) 



[-TT 



Dab 

Qml,AB \ kftT 



(6) 



where k B is the Boltzmann s constant; h denotes the Planck s 

constant; T is the local temperature; yu = is the reduced 

m A B 

mass of the molecule; the 2; nt 's are the internal partition func- 
tions; and Dab is the dissociation energy of the (A, B, AB) reac- 
tion given by Eq. Species A and B themselves can be either 
molecules or single atoms. 

In the identical framework we can consider positive ionisa- 
tion and recombination processes: 

A+ r <_ A+ r+l + £ - 

Again, species A is taken in the general sense and can be either 
a molecule or a single atom, and the superscript +r (or +r + 1) 
denotes its ionisation stage. 

These processes can be described through the corresponding 
equilibrium or ionisation constant: 



n A + 



which is explicitly given in the form of the Saha equation 



2mn„kT 



int,A 



r+ i ■ 2 



h 2 J fiint.A 



exp 



k*T) 



(7) 



(8) 



Here m e is the mass of the electron; I A+r is the ionisation potential 
of species A in the +r th ionisation stage; the Q mt are the internal 
partition functions appropriate to the corresponding species. The 
factor 2 is the statistical weight g e for free electrons, correspond- 
ing to two possible spin states. 

The same formalism with r = — 1 can be applied to account 
for the electron-capture negative ionisation 

A + e~ ±5 A" , 

which is assigned the equilibrium constant 

n A - 



K^ ha (T) = 



and the Saha equation: 



h 2 



2nm e k B T 



3/2 



n A n e - 



(2int,A- / h- 

2int,A-2 \k B T 



(9) 



(10) 



where Ia- corresponds to the electron affinity, i.e. the energy 
released when an electron is attached to a neutral atom or 
molecule. 

Where ionisation of diatomic and polyatomic molecules is 
considered, there are at least three energy-equivalent ways of 
forming positive molecular ions: 



1) A + B-^AB ; 

2) A - e~ -> A + ; 

3) fl - e~ -> B + ; 



AB-e- 
A + + B 
A + B + 



*AB + 
AB + 
AB + . 



Dissociation and ionisation equilibrium can be taken into ac- 
count simultaneously by choosing that dissociation path in 
which the atomic species that remains ionised is the one with 
the lowest ionisation potential. For instance, for a ionised di- 
atomic molecule AB + with I a < Ib the selected sequence is 2), 
so that the number density of the ionised molecule is calculated 
by combining Eq. (O and Eq. ([8J, obtaining: 



Kab+(T) 



v foal 



Saha 



n AB + 



K 



Saha 



AB 



(2nfik B T\ i/2 Qim,A+Qm,B 

bH -q — exp 



int,AB+ 



k B T 



(ID 
(12) 
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where the dissociation energy is given by D A b* - D AB + Ia- Iab 
and Iab is the ionisation energy of the molecule AB. 

In the case of negative molecular ions and assuming that dis- 
sociation of AB produces A~ and B (hence Ig- < Ia), we can 
extend the same formalism of Eq. ( TTZl i to calculate the dissocia- 
tion constant: 



K AB -(T) 



riA-riB 
riAB- 



Saha 



K 



Saha 



AB- 



1 2npk B T\ 3/2 QimA- Qim,B 



Da 



(13) 
(14) 



where the dissociation energy is now Dab- - Dab + Iab- - Ia-, 
and Iab- denotes the electron affinity of AB, or equivalently, the 
neutralisation energy of AB . 



where the summation includes all molecules and atoms (neutral 
and ionised). The number density of each atomic species, N a , is 
then obtained from Eq. (T5[ once the fractions s a = N a /N- d are 
given as a part of the problem specification. 

The foregoing set of equations ((5]) through (fTTT i are sufficient 
for problem solution, as illustrated in the following. 

2.1 .3. Solution to the ICE problem 

The solution to the chemical equilibrium problem in AiSOPUS is 
based in large part on source code available under the GPL from 
the SSynth project (http://sourceforge.net/projects/ssynth/) that 
is developed by Alan W. Irwin and Ana M. Larson. Basic ther- 
modynamic data together with a few FORTRAN routines were 
adopted with the necessary modifications, as detailed below. 



2.1.2. Conservation relations 

In addition to the equilibrium relations (dissociation- 
recombination and ionisation), there exist three additional 
types of equations that will completely determine the con- 
centrations of the various species of the plasma, namely: i) 
conservation of atomic nuclei for each chemical species, ii) 
charge neutrality, and iii) conservation of the total number of 
nuclei. 

Let us denote with A4t the number of chemical elements, 
Mnoi the number of molecules (neutral and ionised), and A/tot the 
total number of species under consideration (neutral and ionised 
atoms and molecules). 

Indicating with N a the number density of nuclei of type a 
(occurring in atoms, ions and molecules), and s a = N a /N a its 
fractional abundance with respect to the total number density of 
nuclei N a (both in atoms and bound into molecules), then the 
conservation of nuclei requires that each atomic species a (not a 
molecule) fulfils the equation 

PZ JVmol 

s a N R — N a — n a + 2^ n a +, + n a - + ^ v A , a n A (15) 

r=l A=\ 

In the right-hand side member, n„ is the number density of neu- 
tral atoms; the next two terms give the number density of ions in 
all positive ionisation stages (up to the maximum stage pz), and 
in the negative ionisation stage; the last summation is performed 
over all molecules (neutral and ionised) which contain the atom 
a. Here va, u corresponds to the stoichiometric coefficient, ex- 
pressing the number of atoms a in molecule A. 
Charge neutrality requires that 

Mot PZ Mot 

n e =2^2^rnA»-^ j n A - (16) 

i=\ r=\ j=\ 

where we include all appropriate atomic and molecular ions, 
with both positive and negative electric charges. For each species 
A,, the total number of free electrons is evaluated by means of 
the second internal summation extended up to pz, which corre- 
sponds to the highest positive ionisation stage. Negative ionisa- 
tion produces a loss of free electrons, which explains the minus 
preceding the last summation. 

Finally, the necessary normalisation is given by the ideal gas 
law, so that the total number density n tot of all particles obeys the 
relation: 

Mot p 

«tot = n e + V n A = — - (17) 



2.1.4. Thermodynamic data 

From the SSynth package we make use, in particular, of the 
whole compilation of internal partition functions, ionisation and 
dissociation energies. Each species (atomic and molecular) is as- 
signed a set of fitting coefficients of the polynomial form 

in 

]nQ = Y j a i QnT) i , (18) 

!=0 

based mostly on the works by Irwin (1981, 1988) and Sauval 
& Tatum (1984). In most cases the degree of the polynomials is 
five (m = 6). The original compilation was partially modified 
and extended to include additional ionisation stages for atoms, 
and two more molecules, Hi and FeH, that may be relevant in 
the opacity computation. We consider the ionisation stages from 
I to V for all elements from C to Ni (up to VI for O and Ne), 
and from I to III for heavier atoms from Cu to U. Specifically, 
our revision/extension of the original Irwin's database involve 
the following species. 

The partition functions for the C to Ni group have been 
re-calculated with the routine pfsaha of the ATLAS 12 code 
(Kurucz 1993a), varying the temperature from 5000 to 20000 K 
in steps of 100 K. The partition functions of the 15 rare earth el- 
ements belonging to the Lanthanoid series, from La to Lu, have 
been re-computed with the routine pfword from the UCLSYN 
spectrum synthesis code (Smith & Dworetsky 1988) increment- 
ing the temperature from 6000 to 30000 K in steps of 100 K. 
This revision was motivated by the substantial changes in the 
energy levels of the earth-rare elements introduced in more re- 
cent years (Alan Irwin, private communication; see e.g. Cowley 
et al. 1994). We have verified that, the UCLSYN partition func- 
tions for third spectra of the Lanthanides are in close agreement 
with the data presented in Cowley et al. 1994, while the results 
from ATLAS12 or Irwin's (1981) compilation are usually lower, 
in some cases by up to a factor of two (e.g. for Ce +3 and Tb +3 ). 
The partition function for FeH is given from Dulick et al. (2003) 
over a temperature range from 1000 to 3500 K in steps of 100 
K. Then, for all the revised species, we have obtained the fitting 
coefficients of Eq. ( fT8l by the method of least-squares fitting. 
In most cases the best fitting is achieved with a x 1 parameter 
lower than 10~ 4 . For H3 we use the original fitting polynomial 
provided by Neale & Tennyson (1995). 

In total, our database of partition functions consists of Wtot ~ 
800 species, including m 300 atoms (neutral and ionised) from 
H to U, and N mo i ~ 500 molecules. 
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Process 


Symbol 


Reaction 


References and Comments 


Rayleigh 


0"Ray(H 2 ) 

0"Ray(H) 

cr Ray (He) 


H 2 + hv — > H 2 + Ay' 
H + hv -> H + Av' 
He + Ay — > He + hv' 


Dalgarno & Williams (1962) 


Dalgarno (1962) 


Thomson 


Th(e-) 


e~ + hv — > e~ + Ay' 


NIST (2006 CODATA recommended value) 


free-free 


o-ff(H-) 
o-ff(H) 

o- ff (H^) 

o- ff (H 3 ) 

o-ff(He-) 

o- ff (He) 

tr ff (He + ) 


H + e~ + Av — > H + e~ 
H + + e" + hv -> H + + e" 

XJ+ , If , 7 . TT+ | TT 

rl + n t IIV — > 1 1 + rl 

H 2 + e~ + hv — > H 2 + e~ 
H, + e~ + hv — > H} + e~ 
He + e~ + hv -» He + e" 
He + + e~ + hv — > He + + e~ 
He ++ + e" + -> He ++ + e" 


John (1988) 

Method as in Kurucz (1970) based on Karsas & Latter (1961) 
Leueaev et ai. (zvuj ) 

Tr»tin HOTM 

0"fj(H3) = <x ff (H) (assumed) 

\~al UOIl CL dl. [ LyQy ) 

o~ ff (He) — o~ ff (H) (assumed) 
trff(He + ) = (Tff(H) (assumed) 


bound-free 


o-bf(H-) 

o-bf(H) 

o-bf(Hj) 

cr bf (He) 

o- b f(He + ) 


H + hv -» H + e" 
H + Ay -> H + + e" 
H} + hv -> H + + H 
He + Ay — » He + + e~ 
He + + Ay -> He ++ + e - 


JOnn \lyyo) 


Method as in Kurucz (1970) based on 

vJlIlgCllCIl \Vy\jy) dllU JVdiSaS OL l^dLlCI ^17U1 ) 


Lebedev et al. (2003) 


A/TptTinH ;iQ in TCiinir"7 ( 1 Q70^ hiiQPfl on 

IVlt, L11UU. CIS ill 1 V Ll 1 \A\*. / . yljs t\JJ ULl&t^U v/1 1 

Gingerigh (1964) and Hunger & Van Blerkom (1967) 


Hunger & Van Blerkom (1967) 


bound-bound 


tr bb (H) 


H + hv -» H* 


Kurucz (1970) including Stark broadening 


Collision 

induced 

absorption 


o-cia(H 2 /H 2 ) 
o-ciA(H 2 /He) 
o-cuCH/He) 


H 2 + H 2 + Av — > H 2 + H 2 
H, + He + Av -> H 2 + He 
H + He + Ay -> H + He 


600 K < T < 7000 K 90 rrrT 1 < v < 70000 cm -1 
Borisow et al. (1997) 


1000 K < T < 7000 K, 25 cmr'y <v< 20088 cnT 1 
J0rgensen et al. (2000) 


1500 K < T < 10000 K, 50 cnr 1 < v < 1 1000 cm" 1 
Borisow et al. (2001) 



Table 1. Scattering and absorption processes involving H and He nuclei, considered in this work. 



2.1.5. Method 



First we need to specify the list of atoms, ions and molecules 
which should be considered, together with the values of gas pres- 
sure P, temperature T and chemical abundances e a = N a IN a . 
Then, the code arranges a system consisting of 7V e i + 2 non- 
linear equations for the number densities of neutral atoms n a , 
the total number density of atoms 7V a , and the electron density 
n e . Once these densities are known, the number densities of any 
other ionised and/or molecular species are calculated by solving 
for their concentrations in Eqs. (|5), 0, ©, ( fTTT i. or ( fl~3T > using 
the equilibrium/ionisation constant appropriate for each atom or 
molecule. Given the non-linearity of the equations, the system is 
conveniently solved by using a standard Newton-Raphson iter- 
ative method (Press et al. 1986). Numerical details are given in 
AppendixlAl 

It is worth remarking that the EOS in itiSOPUS can easily 
deal with any chemical mixture, including peculiar cases such 
as zero-metallicity (Z = 0) or hydrogen-free (X = 0) gas. In 
general, no convergence problem has been encountered within 
the assumed ranges of the state variables. 

In place of the gas pressure P, it is also possible to specify 
the gas density p. In this case a second external iteration cycle 
is switched on according to a root-finding numerical scheme. At 
each i th iteration a new value P, is assigned to the pressure and 
the EOS is solved yielding the corresponding p, = f, (jum u ) / k^T , 
where /i is the mean molecular weight in units of atomic mass 
m u . The process is repeated until the difference | log(p,) - log(p)| 
decreases below a specified tolerance 5 p . In our computations we 
adopt 6 P = 10~ 8 , and convergence is reached typically after 3-4 
iterations. 



Species Source & Reference 





C,N,0 






Ne,Na, Mg 




Atoms 


Al,Si,S 
Ar, Ca, Cr 
Mn,Fe,Ni 


OP: Seaton (2005) for log(T) > 3.6 




CI, NI 


Method as in Kurucz (1970) based on 


Atoms 


OI, Mgl 


Peach (1970) and Henry (1970) 




All, Sil 


for log(T) < 3.6 




HF 


LL: Uttenthaler et al. (2008) 




HC1 


LL: Rothman et al. (2005) 




CH 


LL: J0rgensen (1997) 




C 2 


LL:Querci et al. (1974) 




CN 


LL: J0rgensen (1997) 




CO 


LL: Goorvitch & Chackerian (1994) 




OH 


LL: Schwenke (1997) 




SiO 


LL: Rothman et al. (2005) 




TiO 


LL: Schwenke (1998) 


Molecules 


VO 


LL: Alvarez & Plez ( 1998) 


CrH 


LL: Bauschlicher et al. (2001) 




FeH 


LL: Dulick et al. (2003) 




YO 


LL: Littleton (2007) (2001) 




ZrO 


LL: Plez (2007) 




H 2 


LL: Barber et al. (2006) 




HCN 


LL: Harris et al. (2003) 




c 3 


OS: J0rgensen et al. (1989) 




co 2 


LL: Rothman et al. (1995) 




so 2 


LL: Rothman et al. (2005) 




C,H 2 


OS: J0rgensen (1997) 



Table 2. Data sources for the atomic and molecular monochro- 
matic absorption coefficients. Atomic absorption coefficients 
(including both continuum and discrete opacities) are from the 
Opacity Project (OP) database, while molecular absorption co- 
efficients are extracted from either line lists (LL) or opacity sam- 
pling (OS) data. 
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2.2. Opacity 



with a Doppler width A v given by 



In our computations we consider the following continuum opac- 
ity processes 

- Rayleigh scattering, 

- Thomson scattering, 

- Bound-free absorption due to photoionisation, 

- Free-free absorption, 

- Collision-induced absorption (CIA), 

and line opacity processes 

- Atomic bound-bound absorption, 

- Molecular band absorption. 

Denoting with cry(v) the monochromatic cross section 
(in cm 2 ) of the / absorption process (not scattering), the 
monochromatic true absorption opacity and scattering opacity 
per unit mass (in cm 2 g~') are calculated with 



*f(v) 



*r>) 



P 1 



■ (T 



P 



Xv), 



(19) 
(20) 



where nj is the number density of particles of type j, p is the gas 
density, and (1 - e~ hv ^ kBT ) is a correction factor for stimulated 
emission. 

Tables [T]and|2]detail the whole compilation of the scattering 
and absorption processes considered here. 

The monochromatic opacity cross sections for atoms (ex- 
cept for H and He), taken from the OP database, are interpolated 
in frequency, temperature and electron density, according to the 
formalism described in Seaton et al. (1994) and Seaton (2005). 
They include all radiative continuum and discrete opacity pro- 
cesses. Line broadening is taken into account as the result of 
thermal Doppler effects, radiation damping and pressure effects. 

The monochromatic molecular absorption coefficient caused 
by each of the different species included in our code is taken 
from opacity sampling (OS) files produced for the selected fre- 
quency grid (see Sect. 12.2.21 and Appendix 151. that are in most 
cases calculated directly from the corresponding line lists (see 
Table|2]i. The only exceptions are C2H2 and C3 for which we use 
already existing pre-computed opacity sampling data. 

Where line lists are adopted, the absorption cross section of 
a spectral line, involving the bound-bound transition from state 
m to state «, is evaluated with the relation: 



< S (v) = — 



gf 



-E /k 



m e c Q(T) 



> T (l ~ e - hv o /k » T ) (f>(v) (21) 



with e and m e the charge and mass of the electron, c the speed of 
light, hvo the energy of the corresponding radiation, Q(T) the to- 
tal partition function (being the product, Q lrdns Qi nt , of the trans- 
lational and internal partition functions) of the molecular species 
under consideration, £0 the excitation energy of the lower level 
m of the transition, gf the product of the statistical weight g( m ) 
of the level times the oscillator strength /(,„,„) of the transition. 
The correction for stimulated emission is given by the term in 
brackets. The normalised broadening function, <f>(v), for the line 
profile takes into account the effect of thermal broadening and 
non thermal-contribution of microturbolent velocities, according 
to the equation: 

1 -&f 



A v 



(22) 



A V 

A v = — 



2k n T 



(23) 



where m is the mass of the molecule, and £ is the microturbo- 
lent velocity, which is assigned the value 2.5 km/s. More details 
about the input data and the treatment of molecular line opacities 
can be found in Aringer (2000), Lederer & Aringer (2009), and 
Aringer et al. (2009). 

In summary, to generate the molecular OS files directly from 
the line lists, prior to the execution of iESOPUS, we proceed 
as follows. For each value of a selected set of temperatures, (13 
values in the range 600 K 5 T 5= 10000 K), the monochro- 
matic absorption coefficient of a molecular species at a given 
wavelength point, tr^^v), is obtained by adding up the contribu- 
tions of all the lines in the list with the corresponding broadening 
functions taken into account: 



lines 



(24) 



where each term 0"j^ s (v) is evaluated with Eqs. ( |2"TT i - d23l ). Then, 
during the computations with iESOPUS, we interpolate on the 
OS tables for any given temperature of the gas. We notice that 
the errors brought about by this interpolation are marginal com- 
pared to all other sources of uncertainty (e.g. molecular data, 
microturbolence velocity, solar abundances, etc.). 



2.2.1. The Rosseland mean 

Once the total monochromatic opacity coefficient is obtained by 
summing up all the contributions of true absorption and scatter- 
ing 

<v) = J\f(v) + «f tt (v) , (25) 
j 

then the Rosseland mean opacity, classically defined by Eq. (f2j) 
is conveniently calculated with (see e.g. Seaton et al. 1994): 



where 



F R (u) 



kr(P 



15 
(4tt 4 ) 



1 H Fr{u) a 

0, T) J k(u) 



u exp(-w)/[l - exp(-w)] . 



(27) 



In the above equations v is the photon frequency, and 
u = liv/ik-sT) is the normalised photon energy. In our calcu- 
lations k r denotes the absorption coefficient per unit mass, and 
is always given in cm 2 g _1 . Since the opacity coefficient k v en- 
ters Eq. (fJJ as an inverse, the minima dominate the values of 
the Rosseland mean. It follows that a large k r implies large ab- 
sorption from the radiation beam, while a small a-r indicates that 
the energy losses from the beam remain little as it propagates 
through the matter. 

In practise, the numerical integration of Eq. d26b requires to 
specify two finite (lower and upper) limits, u\ and ui, and the 
grid of frequency points. The choice of the limits must guarantee 
the covering of the relevant wavelength region for the weighting 
function dB v /dT, so as to include its maximum and the declining 
wings. 

In this respect it useful to recall that, in analogy with the 
Wien's displacement law for the Planck function, the wavelength 
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A max of the the maximum of dB v /dT is inversely proportional to 
the temperature according to 

3756.56 

Amax [jum] = (28) 
1 Lis. J 

It follows that the maximum of the function Fr is reached for 
Mmax = 3.8300. 

In our calculations we adopt the integration limits U\ - 1 0~ 3 
and M2 — 64, corresponding to the wave numbers v\ = 10 
cm -1 and V2 = 2 10 5 cm -1 , and wavelengths A\ = 1000 
/mi and A2 = 0.04 ^m, respectively. We have verified that 
these values largely satisfy the condition of spectral coverage 
of the weighting function over the entire temperature range, 
3.2 < log(r) < 4.5, here considered. 

2.2.2. The frequency grid and computing time 

Since in our calculations a number of crucial opacity sources, 
i.e. molecular absorption bands, are included as OS data, it is 
convenient to specify, prior of computations, a grid of frequency 
points, which should be common to both the OS treatment and 
the numerical integration of Eq. ©. The frequency distribution 
will be determined as a compromise between the precision (and 
accuracy) of the integration and the speed of calculations. 

For this purpose we employ the algorithm by Helling & 
J0rgensen (1998), that was developed to optimise the frequency 
distribution in the opacity sampling technique when dealing 
with a small number of frequency points. We performed a few 
tests adopting frequency grids with decreasing size, namely with 
n tot = 5488, 1799, 944, 510, and 149 frequency points. More 
details are given in Sect. [B] The results discussed in the follow- 
ing sections refer to the grid with n tot = 944 points, which has 
proved to yield reasonably accurate RM opacities. 

Besides the quality of the results, another relevant aspect is 
the computing time. With the present choice of the frequency 
grid , i.e. n lot = 944 points, generating one table at fixed chem- 
ical composition, arranged with the default grid of the state pa- 
rameters (T and R, see Sect. 13.11 ). i.e. containing NtXNr = 67 X 
19 = 1273 opacity values, takes r ~ 45 s with a 2.0 GHz proces- 
sor. Adopting other frequency grids would require shorter/longer 
computing times, roughly t ~ 200 s for n lot = 5488; t ~ 70 s 
for n tot = 1799; r ~ 30 s for n tot = 510; and t ~ 15 s for 
«tot = 149. These values prove that jESOPUS is indeed a quick 
computational tool, which has made it feasible, for the first time, 
the setup of a web-interface (http://stev.oapd.inaf.it/aesopus) to 
produce low-temperature RM opacity tables on demand and in 
short times. 

The main reason of such a fast performance mainly resides 
in the optimised use of the opacity sampling method to describe 
molecular line absorption, and the adoption of pre-tabulated ab- 
sorption cross-sections for metals (available from the Opacity 
Project website). In this way the line-opacity data is extracted 
(e.g. from line lists and the OP database) and stored in a conve- 
nient format before the execution of iESOPUS, thus avoiding to 
deal with huge line lists during the opacity computations. This 
latter approach is potentially more accurate, but extremely time- 
consuming (e.g. F05). 

Moreover the improvement in accuracy that would be 
achievable with the on-the-fly treatment of the line lists is in 
principle reduced when adopting a frequency grid for integra- 
tion which is much sparser (e.g. ~ 10 4 frequency points as in 
F05) than the dimension of the line lists (up to 10 7 - 10 8 line 
transitions). On the other hand, as shown by our previous tests 




-10 -5 5 10 



log(R) 

Fig. 1. Location of our RM opacity tables in the log(T) - \og{R) 
diagram (shaded rectangular area), together with the approxi- 
mate boundaries between regions where the total pressure is 
dominated by different effects: radiation pressure, ideal gas, 
electron-degenerate gas, Coulomb interactions, and pressure 
ionisation. The vertical line to the left of which radiation dom- 
inates the pressure is given by P gas = P ra( j with fi = 0.5. Along 
the border line between the ideal and degenerate gas we equate 
the corresponding pressures of a non-relativistic electron gas 
with mean molecular weight [i e = 2. The vertical boundary at 
\og(R) ~ 3, beyond which the Coulomb coupling of charged 
particles should become important, is defined by the condition 
Fc = 1, where Tc = 1.1 10~ 5 77p 1/3 is the Coulomb coupling 
parameter for an ionised-hydrogen plasma. Pressure ionisation 
is assumed to become dominant at log(p) = for log(T) < 4.5, a 
typical value according to the analysis developed by Luo (1997). 
The behaviour of the parameter ft (defined in the text) is also 
shown. As an example, the evolution of the photospheric param- 
eters (r e ff, Rphot) of a stellar model with initial mass M = 5 M Q , 
and metallicity Z = Z re f = 0.02 is depicted by a magenta line, 
covering the evolution from the pre-main sequence to the first 
pulses on the TP-AGB (calculations performed with the Padova 
stellar evolution code). 



and also by F05, while the computing time scales almost lin- 
early with the number of frequency points, the gain in precision 
does not, so that the RM opacities are found to vary just neg- 
ligibly beyond a certain threshold (see also Helling et al. 1998 
and Appendix|B]l. All these arguments and the results discussed 
in Sects. 14.1.11 support the indication that the agile approach 
adopted in jESOPUS is suitable to produce RM opacities with 
a very favourable accuracy/computing-time ratio. 



3. Opacity tables: basic parameters 

Tables of RM opacities can be generated once a few input pa- 
rameters are specified, namely: the chemical composition of the 
gas, and the bi-dimensional space over which one pair of inde- 
pendent state variables is made vary. 
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3.1. State variables 

Under the assumption of ideal gas, described by the law 
^gas = pT , 



(29) 



one must specify one pair of independent state variables. Usual 
choices are, for instance, CP gas , T) or (p, T). For practical and 
historical reasons, opacity tables are generally built as a function 
of the logarithm of the temperature T, and the logarithm of the 
R variable, defined as R = p/(T 6 ) 3 , with T 6 = 77(10 6 K). 

An advantage of using the R parameter, instead of p or P, 
is that the opacity tables can cover rectangular regions of the 
(R, T)-plane, without the nasty voids over extended temperature 
ranges that would come out if intervening changes in the EOS 
are not taken into account (e.g. transition from ideal to degener- 
ate gas). 

Interestingly, as pointed out by Mayer & Duschl (2005; see 
their appendix D), different R values correspond to different 
gas/radiation pressure ratios, 6 = P ga s/-P ra d- The relation be- 
tween \og{R) and log(6) is linear, with larger R values corre- 
sponding to larger 6, i.e. an increasing importance of P gas against 
f la d- Moreover, we notice that the equality P gas = P ra d takes 
place in the range at -4.8 ^ log(R) $ -4.5, assuming a mean 
molecular weight varying in the interval 0.5 $ fx £ 1. In Fig.Q] 
we also plot the quantity /3 = P gas /(Pgns + Prad), a parameter 
frequently used by stellar evolutionists. 

In this respect Fig. [TJ illustrates the rectangular region cov- 
ered by our RM opacity tables in the log(T) - \og{R) dia- 
gram, defined by the intervals (3.2 < log(T) < 4.5) and 
(-8 < log(fl) < 1). We note that the table area lies in the 
domain of the ideal gas, and it extends into the region dominated 
by radiation pressure for \og(R) < -4.5. Non ideal effects related 
to electron degeneracy, Coulomb coupling of charged particles, 
and pressure ionisation of atoms are expected to become domi- 
nant outside the table boundaries, in the domain of high-density 
plasmas. 

It is important to remark that our RM opacity tables can be 
easily extended to higher temperatures, log(T) > 4.5, with the 
RM opacity data provided by OPAL and OP. As a matter of fact 
the agreement between our results and OPAL is good in the over- 
lapping transition region, say 3.9 $ log(T) 5 4.5 (see Sect. l4.1.T1 
and panel c) of Fig [7]). 

Within the aforementioned limits of the state variables, the 
interactive web mask enables the user to freely specify the ef- 
fective ranges of log(T) and \og(R) of interest as well as the 
spacing of the grid points A log(T) and A log(R). From our tests 
it turns out that a good sampling of the main opacity features 
can be achieved with Alog(T) = 0.05 for log(T) > 3.7 and 
Alog(r) = 0.01 for log(T) < 3.7, and Alog(fl) = 0.5. In any 
case, the choice should be driven by consideration of two as- 
pects, i.e. maximum memory allocation, and accuracy of the 
adopted interpolation scheme. 

3.2. Chemical composition 

It is specified in terms of the following quantities: 

- The reference solar mixture; 

- The reference metallicity Z re f ; 

- The hydrogen abundance X; 

- The reference mixture; 

- The enhancement/depression factor f, of each element (heav- 
ier than helium), with respect to its reference abundance. 



Reference 




(C/(J)o 


(f~* tr\\Q 
(C/U£prit>l 


Anders & Orevesse l^oV (ALroy) 


A A 1 A A 


0.42/ 


a nco 






0.479 




Grevesse & Sauval 1998 (GS98) 


0.0170 


0.490 


0.947 


Holweger 2001 (H0l£, 


0.0149 


0.718 


0.937 


Lodders 2003 (L03) 


0.0132 


0.501 


0.929 


Grevesse et al. 2007 (GAS07) 


0.0122 


0.537 


0.929 


Caffau et al. 2009 (C09£ 


0.0155 


0.575 


0.938 



" This abundance ratio is defined by Eq. d35t . 

* The elemental abundances are taken from Grevesse & Sauval (1998), 
but for C, N, O, Ne, Mg, Si, and Fe that are modified following the revi- 
sion by Howeger (2001). 

c The elemental abundances are taken from Grevesse & Sauval (1998), 
but for N, O, and Ne following the revision by Caffau et al. (2008, 2009). 

Table 3. Compilations of the solar chemical composition 
adopted in the computation of the EOS and gas opacities. For 
each mixture the solar total metallicity Z Q (in mass fraction), the 
abundance ratios (C/O) and (C/0) cr jt,i are indicated for compar- 
ison. The latter marks a critical boundary for the gas molecular 
chemistry in the range 3.2 < log(T) < 3.6. The C and O 
abundances are expressed as number fractions. 



The reference solar mixture can be chosen among various 
options, which are referenced in Table [3] For their relevance to 
the opacity issue, the corresponding solar metallicity, Z Q , and 
the (C/O) raticQ are also indicated. Scrolling Table[5]from top 
to bottom we note that Z significantly decreases, passing from 
~ 0.019 in AG89 down to ~ 0.012 in GAS07. This implies that 
opacity tables constructed assuming the same Z may notably dif- 
fer depending on the adopted solar mixture. Concerning C/O, a 
key parameter affecting the opacities for log(T) $ 3.5, we see 
that it spans a rather narrow range (0.43 ^ C/O ^ 0.53) pass- 
ing from one compilation to the other, except for the H01 which 
corresponds to a higher value, C/O ~ 0.72. How much these 
differences in the reference solar mixtures may impact on the 
resulting opacities is discussed in Sect. 14.11 

Let us indicate with Nz the number of metals, i.e. the chem- 
ical elements heavier than helium, with atomic number Z, ^ 3. 
Each metal is characterised by an abundance X, in mass fraction 
and, equivalently, an abundance e, in number fraction, respec- 
tively defined as: 



X, 



AiNt 



and E: = — 



(30) 



where 7Y, is the number density of nuclei of type i with atomic 
mass A,, and N a is the total number density of all atomic species 
(with the same notation as in Sect. 12.1.2b . In both cases the nor- 
malisation condition must hold, i.e. ^.Jj Xj = 1 and Y!^t\ s i - !• 
The total metal abundance is given by Z = 2;>3 Xj in mass frac- 
tion, and Ez = 2/>3 s i m number fraction. 

We assign each metal species the variation factors, f, and gj, 
relative to the reference mixture: 

Xi = fiX ijet and s t = g,e,\ re f . (31) 

The reciprocal relations between X/ and e, derive straightfor- 
wardly: 



Hp?] A j E j 



and e, = 



(32) 



1 Throughout the paper the C/O ratio is calculated using the abun- 
dances of carbon and oxygen expressed as number fractions, i.e. C/O = 
ec/so following the definition given by Eq. < |301 >. 
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Fig. 2. Fractional abundances of elements, with nuclear charge 
Zj = 6 — 30, normalised to the solar metallicity according to 
various compilations, as indicated. 



log(K, 




3.6 3.8 
log(T) 



Fig. 3. Rosseland mean opacity as a function of variables T and 
R over the entire parameter space considered in our calculations. 
The adopted composition is assumed to have Z = Z, e f = 0.02, 
X = 0.7 and the metal abundances scaled-solar to the GS98 mix- 
ture. 



as well as those between fi and gj for metals: 



fi = gi z and 

2j j>3 A j s j 



8i = fi 



^>3 



X 



Zj>3Xj/Aj 



(33) 



We have verified that /, ss g, as long as they are not too large 
and the ratios between the two summations in the left-hand side 
members of Eq. ( T33l > do not deviate significantly from unity (see, 
for instance, Table 0). 

In principle, the reference chemical mixture can be any given 
chemical composition. Frequent choices are, for instance, mix- 
tures with scaled-solar partitions of metals, or with enhanced 
abundances of a-elements. The jESOPUS code is structured to 
allow large freedom in specifying the reference mixture. For 
simplicity, in the following we will adopt the solar mixture as the 
reference composition, so that the reference metal abundances 



and 



(34) 



with clear meaning of the symbols. The partitions, X, i0 /Z , of 
chemical elements from C to Zn are shown in Fig. [2] for a few 
compilations of the solar chemical composition. 

According to the notation presented by Annibali et al. 
(2007), the chemical elements can be conveniently divided into 
three classes depending on the sign of /; (or gj) , namely: 

- enhanced elements with jj > 1 (org,- > 1); 

- depressed elements with fi < 1 (org; < 1); 

- fixed elements with fi = 1 (org/ = 1). 

The latter correspond to the reference abundances, i.e. scaled- 
solar in the case discussed here. Moreover, let us designate 

- selected elements with fi + 1 (org,- + V) 

the group of elements which are assigned variation factors dif- 
ferent from unity (either enhanced or depressed), as part of the 
input specification. We limit the discussion here to the case of 
the abundances X t expressed in mass fraction, since exactly the 
same scheme, with the due substitutions, can be applied to the 
abundances e, in number fraction. In this respect one should bear 
in mind that the conversions X, ^± e, are obtained with Eqs. d32l . 
Starting from the reference mixture, then the new mixture can be 
obtained in two distinct ways: 



1. 



Case Z + Zref . The enhancement/depression factors fi of the 
selected elements produce a net increase/depletion of total 
metal content relative to the reference metallicity Z re f. The 
actual metallicity is calculated directly with Z = Yr=\ fi^itef- 
In this case all Nz variation factors fi can be freely specified 
without any additional constrain. 

Case Z = Zref. The enhancement/depression factors produce 
non-scaled-solar partitions of metals, while the total refer- 
ence metallicity Z re f is to be preserved. This constraint can 
be fulfilled with various schemes, e.g. by properly varying 
the total abundance of all other non-selected elements so as 
to balance the abundance variation of the selected group. For 
instance, if the selected elements have all fi > 0, so that 
we refer to them as enhanced group, then the whole posi- 
tive abundance variation should be compensated by the neg- 
ative abundance variation of the complementary depressed 
group. Another possibility is to define a fixed group of el- 
ements whose abundances should not be varied, hence not 
involved in the balance procedure; in this case the preserva- 
tion of the metallicity is obtained by properly changing the 
abundances of a lower number of atomic species among the 
non-selected ones. 

In principle, the quantities fi can be chosen independently for 
up to a maximum of (Nz - 1) elements, while the remaining 
factor is bound by the Z = Z re f condition. A simple practise is 
to assign the same factor to all the elements belonging to the 
selected group, either enhanced or depressed, as frequently 
done for c-enhanced mixtures. In this respect more details 
can be found in Sect. 14.31 



The former case (Z + Z e f) properly describes a chemical 
mixture in which the abundance variations are the product of 
nuclear burnings occurring in the stellar interiors. This applies, 
for instance, to thermally-pulsing asymptotic giant branch (TP- 
AGB) stars whose envelope chemical composition is enriched 
in C and O produced by He-shell flashes and convected to the 
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surface by the third dredge-up, which leads to an effective incre- 
ment of the global metallicity (Z > Z re f). 

The latter case (Z = Z re f) corresponds, for instance, to 
chemical mixtures with a scaled-solar abundance of CNO ele- 
ments XcnO' but different ratios e.g. Xc/Xqno, Xj^/Xcno, and 
Xo/Xcno- Alternatively, if we consider the abundances in num- 
ber fractions, the condition, scno = const., may describe the sur- 
face composition of an intermediate-mass star after the second 
dredge-up on the early AGB, when products of complete CNO- 
cycle are brought up to the surface. In this case the total number 
of CNO catalysts does not change, while C and O have been 
partly converted to 14 N. Another example may refer to a- en- 
hanced mixtures with different [ff/Fe]> but the same metal 
content Z. 

Finally, it should be noticed that, once the actual metal- 
licity Z is determined, in both cases the normalisation condi- 
tion implies that the helium abundance is given by the relation 
Y = 1 - X - Z. 

4. Results 

In the following sections we will discuss a few applications of 
the new opacity calculations, selecting those ones that may be 
particularly relevant in the computation of stellar models. For 
completeness, our results are compared with other opacity data 
available in the literature. 

4.1. Scaled-solar mixtures 

Let us first illustrate the case of scaled-solar mixtures, which 
will serve as reference for other compositions. As an example, 
Fig. [3] visualises the tri-dimensional plot of one opacity table 
calculated over the whole log(T) - log(R) parameter space for a 
given chemical mixture. The latter is characterised by (X = 0.7; 
Zref = 0.02; Z = Zref; fi = 1, for i = 3 ■ ■ ■ N e \) according to the 
notation introduced in Sect. [3] meaning that all metal abundances 
are scaled-solar. One can see that the grid of the state variables 
(i.e. AlogfJ) = 0.01 for 3.2 < log(T) < 3.5, and Alog(r) = 
0.05 for 3.5 < log(7/) < 4.5; Alog(#) = 0.5) is sufficiently dense 
to allow a smooth variation of kr all over the parameters space, 
which is a basic requirement for accurate interpolation. 

Different opacity sources dominate the total kr in different 
regions of the log(T) - \og(R) plane. Roughly speaking, we may 
say that the continuous and atomic opacities prevail at higher 
temperatures, while molecular absorption plays the major role 
for log(T) 55 3.5. It has been known for long time (see e.g. 
Alexander 1975), for instance, that the prominent opacity bump 
peaking at log(T) =~ 3.25 in Fig.[6]is mainly due to the strong 
absorption of H2O molecular bands. To delve deeper into the 
matter it is instructive to look at Fig. [4] and Fig. [5] which il- 
lustrate the basic ingredients affecting the RM opacity and their 
dependence on wavelength, temperature and density. 

Figure [4] displays the spectral behaviour of the monochro- 
matic opacity coefficient per unit mass, K/(v), of several ab- 
sorption and scattering processes, as defined by Eqs. (fT9b — (f20b- 
We consider three representative values of the temperature (i.e. 
log(T) = 3.3, 3.7, 4.0) and three choices of the R variable (i.e. 
\og(R) = -8, -3, 1), for a total of nine panels that should 
sample the main opacity domains. For each temperature, we 
also indicate in Fig. [4] the spectral range most relevant for the 
Rosseland mean, by marking the wavelength, A max , at which the 
Rosseland weighting function reaches its maximum value (given 
by Eq.l28li. and the interval across which it decreases by a factor 
He. 



At larger temperatures, i.e. \og(T) = 4.0 and /} max ~ 0.38 pm 
(top panels), the total monochromatic coefficient is essentially 
determined by the Thomson e~ scattering at very low gas densi- 
ties (see the top-left panel for \og(R) = -8), while the H opacity 
(bound-bound, bound-free, and free-free transitions) plays the 
major role at large p. Next to hydrogen, some non-negligible 
contribution comes from atomic absorption at shorter wave- 
lengths. 

At intermediate temperatures, i.e. log(T) = 3.7 and A m . dx ~ 
0.75 pm (middle panels), Thomson e~ scattering again controls 
the total absorption coefficient at the lowest densities, whereas 
at increasing p the most significant opacity sources are due to 
metals and H absorption (electron photo-detachment for A < 
1 .644 pm and free-free transitions). 

At lower temperatures, i.e. log(T) = 3.3 and A max ~ 1.88 pm 
(bottom panels), the molecular absorption bands (mainly of 
H2O, VO, TiO, ZrO, CO) dominate the total absorption coef- 
ficient at any gas density except for very low values, where the 
spectral gaps between the molecular bands are filled in with the 
Thomson e~ scattering coefficient. Due to its harmonic charac- 
ter, the Rosseland mean opacity emphasises just these opacity 
holes, so that the total kr for log(r) = 3.3 and \og(R) = -8 
will be mostly determined by the Thomson e~ scattering, with a 
smaller contributions from molecules. 

This fact becomes more evident with the help of Fig. [5] 
which provides complementary information on both the chem- 
istry of the gas, and the characteristic temperature windows of 
different opacity sources. Results are presented as a function of 
temperature for three values of the parameter R. 

As for the chemistry (top panels of Fig. [5j, we show the con- 
centrations of a few species, selecting them among those that 
are opacity contributors, while leaving out all other chemicals 
to avoid over-crowding in the plots (we recall that iFiSOPUS 
solves the chemistry for Ni ot ~ 800 species). It is useful to re- 
mark a few important features, namely: i) at lower temperatures 
molecular formation becomes more efficient at increasing den- 
sity, ii) the most abundant molecule is either carbon monoxide 
(CO) thanks to its high binding energy at low and intermediate 
densities, or molecular hydrogen (H2) at higher densities; iii) the 
electron density n e is essentially supplied by H ionisation down 
to temperatures log(T) ^ 3.8 - 3.6, below which the main elec- 
trons donors are nuclei with low-ionisation potentials, such as: 
Mg, Al, Na, Si, Fe, etc. (see Fig.|22]and Sect.|43]for more dis- 
cussion of this point). 

The bottom panels display the contributions of several ab- 
sorption/scattering processes to the total RM opacity. This is 
done by considering, for a given source j, the ratio /e° ff //cr, where 

k oS is the reduced RM opacity obtained by including all opacity 

sources but for the / h itself. 

At very low densities, i.e. \og(R) = -8 (left-hand side panel 
of Fig. |5]l the most important opacity source, all over the tem- 
perature range under consideration, is by far Thomson scattering 
from free electrons. Note that at lower temperatures a relatively 
important contribution is provided by Rayleigh scattering from 
neutral hydrogen, while the role of molecules is marginal since 
at these low densities molecular formation is inefficient. 

Different is the case with \og(R) = -3 (middle panel of 
Fig. [5]). We can distinguish three main opacity domains as a 
function of temperature. At lower temperatures, say for 3.2 % 
log(T) ^ 3.6, molecules completely rule the opacity, with H2O 
being the dominant source for log(T) 5 3.4. Additional mod- 
est contributions come from metal oxides, such as TiO, VO, 
YO, and SiO. Note that, though for C/O < 1 the chemistry 
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Fig. 4. Monochromatic absorption coefficients for several opacity sources as a function of the wavelength, for three values of the 
temperature and three values of the R variable, as indicated. The chemical composition is defined by Z = Z le f = 0.02, X = 0.7 
with metal abundances scaled-solar to the GAS07 mixture. The total coefficient is depicted by the highest black line. The vertical 
arrow marks the wavelength of the maximum of the Rosseland weighting function (given by Eql28T>. while the horizontal arrows 
delimit the wavelength range within which the Rosseland weighting function drops by a factor 1/e. Where molecular absorption 
bands are important, the corresponding spectral intervals are also indicated. For graphical purpose only, line absorption coefficients 
for molecules and atoms are smoothed by convolution with a Gaussian function. The variance is empirically chosen to depend on 
the wavelength so as to have a neat representation without missing important spectral details. 



is dominated by O-bearing molecules, there is a small opacity 
bump due to CN at log(r) as 3.5. At intermediate temperatures, 
3.6 ^ log(T) 3.8, the most important role is played by the 
H~ continuum opacity, which in turn depends on the availability 
of free electrons supplied by ionised metals. Additional opac- 
ity contributions are provided by Thomson scattering from elec- 



trons and Rayleigh scattering from neutral hydrogen. At larger 
temperatures, 3.8 £ log(r) ^ 4.5, the total RM opacity is de- 
termined mostly by the b-f and f-f continuous absorption from 
hydrogen, with further contributions from b-b transitions of H 
and atomic opacities. 
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Fig. 5. Top panels: concentrations of various chemical species as a function of temperature, for three values of the R parameter, as 
indicated. Bottom panels: contributions of different opacity sources (both continuous and line-absorption processes) to the total RM 
opacity. Each curve corresponds to log(/c)-log(K° ff ), where k is the full opacity including all opacity sources here considered, whereas 

K ofl is the reduced opacity obtained omitting one particular source at once (labelled nearby). The logarithmic notation allows to 
highlight the temperature domains which characterise the different opacity contributors. The adopted chemical composition consists 
of X = 0.7, Z = Z re f = 0.02, with elemental abundances scaled according to the GAS07 solar mixture. 



In the high density case with log(R) = 1 (right-hand side 
panel of Fig. [5j, the opacity pattern is similar to the one just 
described, with a few differences. The most noticeable ones are 
the sizable growth of the H opacity bump in the intermediate 
temperature window, and the increased importance of the H lines 
at higher temperatures. 

Finally, we close this section by examining the sensitiveness 
of the RM opacity to the underlying reference solar mixture. 
Figure [6] shows an example of our opacity calculations made 
adopting a few solar abundances compilations available in the 
literature. They are summarised in Table [3] The largest differ- 
ences are expected for logT S5 3.4, where the RM opacity is 
dominated by the opacity bump caused by the H2O molecule, 
whose amplitude is extremely sensitive to the excess of oxygen 
with respect to carbon, hence to the C/O ratio. In fact, we no- 
tice that the opacity curves corresponding to GN93, GS98, L03, 
GAS07, and C09 lie rather close one to each other, just reflecting 
the proximity of their C/O ratios (=s 0.5 - 0.6; see Table [3). For 
the same reason, the RM opacity predicted at log T ~ 3.3 with 
the H01 solar mixture is roughly 50% lower, given the higher 
C/O ratio (« 0.7). 

Some differences in RM opacity are also expected in the 
3.5 ^ log(T) ^ 3.65 interval, which is affected mainly by the CN 



molecular bands and the negative hydrogen ion H . We see in 
Fig.|6]that most of the results split into two curves: the opacities 
based on L03 and GAS07 (and partly also C09) are higher than 
those referring to GN93 and GS98 solar mixtures. In this case 
the differences are not caused by the CN molecule, but rather 
reflect the differences in the electron density. As one can notice 
in Fig. |2] L03, GAS07 (and C09) compilations correspond to 
higher solar partitions, X, G /Z , of those elemental species that 
mostly provide the budget of free electrons at these tempera- 
tures, such as: Mg, Si, Ca, and Fe (see also Fig. [22b . As a con- 
sequence, the HF opacity is strengthened in comparison to the 
GN93 and GS98 cases. On the other hand, the opacity curve 
corresponding to the H01 mixture lies somewhere in the middle. 
This is the indirect result of the larger C/O ratio (i.e. more carbon 
is available) which favours a larger concentration, hence opacity 
contribution, of the CN molecule in this temperature window. 



The arguments developed here indicate that the expression 
"standard solar composition" should be always specified explic- 
itly together with its reference compilation and not taken for 
granted, since significant differences arise in the RM opacities 
depending on the adopted solar mixture. 
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Fig. 6. Rosseland mean opacity as a function of temperature 
and assuming log/? = -3. The adopted chemical composition 
consists of X = 0.7, Z = Z re f = 0.02, with elemental abun- 
dances scaled according to a few compilations of the solar mix- 
ture abundances, namely: Grevesse & Noels 1993; Grevesse & 
Sauval 1998; Holweger 2001; Lodders 2003; Grevesse, Asplund 
& Sauval 2007; Caffau et al. (2008, 2009). Note the significant 
depression of the H2O bump in the Holweger (2001) case com- 
pared to the others, due to the lower oxygen abundance, hence 
(C/O) ratio. 



4.1 .1 . Comparison with other authors 

As a next step we checked our opacity results against tabu- 
lated RM data made publicly available from other authors. In 
Fig.|7]we show eight representative comparisons, based on: the 
widely-used and well-tested database set up by the Wichita State 
University group, i.e. Alexander & Ferguson (1994), Ferguson 
et al. 2005 (hereafter also F05); the recent data by Lederer & 
Aringer 2009 (hereafter also LA09) stored in the VizieR ser- 
vice; the RM data available in the Robert L. Kurucz' home- 
page, and the OPAL and OP data computed via their interactive 
web-masks. The R and T intervals are different depending on the 
source considered. For instance, the comparisons with the OPAL 
and OP opacities cover the range from 3.8 < log(T) < 4.5, since 
no molecular contribution is included in the OPAL and OP data. 

In general we can conclude that the check is quite satisfac- 
tory in all cases under examination, as our opacity values agree 
with the reference data mostly within +0.05 dex, with the largest 
differences reaching up to w ±0.10-0.20 only in narrow regions. 

Let us start discussing the comparison with Alexander & 
Ferguson (1994) and Ferguson et al. (2005), illustrated in pan- 
els from a) to d) assuming various reference solar compositions. 
First we notice that the small magenta areas in the upper-left 
corners of the four panels are not included in the test, since at 



those densities and temperatures dust is expected to condensate 
0, whereas our EOS describes the matter in the gas phase. 

Besides this, in all cases the agreement between the opacity 
data of the Wichita State University group and jESOPUS is very 
good for 3.4 < log(J) < 4.5, the differences Alog(/ffi) being 
mostly comprised within ±0.05 dex throughout the R range. For 
log(T) < 3.4 the deviations between F05 and iESOPUS appear 
to grow with a systematic trend, i.e. log(/<^ sopus ) > log(^ 05 ), 
at increasing R. Anyhow, the variations are not dramatic, the 
biggest values arriving at « -0.15/ - 0.20. This result is not 
surprising since this is just the region where molecular absorp- 
tion dominates, so that the predicted RM opacity is sensitive to 
differences in the treatment of the molecular line opacities (line 
lists, broadening, adopted frequency grid, etc.). 

This applies also when comparing different releases of the 
same database as it is illustrated, for instance, by panels a) and 
b) relative to the data of the Wichita State University group. We 
notice that where jESOPUS exhibits the best agreement (< 0.05 
dex) with Alexander & Ferguson (1994) at log(T) * 3.4 and 
\og(R) ^ -3, the largest differences (0.15 - 0.20 dex) show up 
instead in the comparison with F05 for the same set of abun- 
dances. In this respect, we expect that much of the discrepancy 
between F05 and jESOPUS for 3.2 < log(r) < 3.4 is due to 
the different molecular line data adopted for water vapour, i.e. 
Partridge & Schwenke (1997) and Barber et al. (2006), respec- 
tively. 

Support to the above interpretation is found when compar- 
ing panel c) and e), the latter showing the check of iESOPUS 
results against Lederer & Aringer (2009) for the L03 solar mix- 
ture. As we see the agreement here is quite fair all over the 
log(T) - log(R) diagram, even in the low-T corner dominated 
by H2O, VO, and TiO absorption, where larger differences with 
F05 (panel c) arise. As a matter of fact, in iESOPUS we adopt 
essentially the same molecular data as in LA09, so that a good 
match is in principle expected. 

Finally, let us briefly comment on the bottom panels (g and 
h) of Fig. [7] relative to two data sets, OP and OPAL, which are 
widely used to describe the RM opacity of the gas in the high-T 
regions, say for log(T) > 4.0. The comparison with iESOPUS 
in the overlapping interval, 3.8 < log(T) < 4.5, is really excel- 
lent, so that the OP and OPAL opacity tables may be smoothly 
complemented in the low-T regime with iESOPUS calculations. 

4.1.2. Tests with stellar models 

The numerical differences in a-r between different authors, illus- 
trated in previous Sect. I4.1.1I assume a physical meaning when 
one analyses their impact on the models in which the Rosseland 
mean opacities are employed. As already mentioned in Sect. [1] 
the largest astrophysical use of pre-tabulated kr(p, T) is in the 
field of stellar evolution models to describe, in particular, the 
thermodynamic structure of the most external layers including 
the atmosphere. 

While it is beyond the scope of this paper to perform a de- 
tailed analysis of the effects of low-T opacities on stellar struc- 
ture and evolution, we consider here two illustrative cases, i.e. 
the predicted location in the H-R diagram of the Hayashi tracks 
described by low-mass stellar models while evolving through 
the RGB and AGB phases. To investigate the differences in r e ff 
brought about by different choices of low-T opacity tables, we 



2 The inclusion of dust in pre-computed opacities is in any case prob- 
lematic since in real stars it will hardly form under equilibrium condi- 
tions. 
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Fig. 7. Comparison between our RM opacity results and those provided by other authors, in terms of log(/f™ thor ) - log(/<-^ sopus ). 
Contour lines, with an incremental step of 0.05 dex, are superimposed to guide the eye. In all cases, except for Kurucz , the adopted 
chemical composition corresponds to Z = Zref = 0.02, X = 0.7. External data are taken from: Alexander & Ferguson (1994) and 
Ferguson et al. (2005) adopting the Grevesse & Noels (1993) solar mixture (panels a and b); Ferguson et al. (2005) assuming the 
solar abundances from Lodders (2003) (panel c) and Grevesse & Sauval (1998) (panel d); Lederer & Aringer (2009) adopting the 
Lodders (2003) solar mixture (panel e); Kurucz' web database for a chemical composition with Z = Z G = 0.0194, X = X = 0.7065 
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Fig. 8. Predicted RGB tracks described by a 1 .5 M model with Z = 0.02, X — 0.7 and scaled-solar abundances of metals according 
to either GN93 (left panel) or L03 (right panel). The luminosity is derived from the core-mass luminosity relation given by Boothroyd 
& Sackmann (1988), while increasing the core mass from 0.20 M to 0.45 M Q . The effective temperature is the result of envelope 
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Fig. 9. Predicted AGB tracks described by a 2.0 M model with Z = 0.02, X = 0.7 and scaled-solar abundances of metals according 
to either GN93 (left panel) or L03 (right panel). The luminosity is derived from the core-mass luminosity relation, at the quiescent 
stage of the pre-flash maximum, given by Wagenhuber & Groenewegen (1998), while increasing the core mass from 0.50 M 
to 0.75 M . The effective temperature is the result of envelope integrations (see the text for more details). The different curves 
correspond to different RM opacity tables, in the temperature range 3.2 < log(r) < 4.0. 



have carried out numerical integrations of a complete envelope 
model (basically the same as the one included in the Padova stel- 
lar evolution code) which extends from the atmosphere down to 
surface of the degenerate core. The overall numerical procedure 
is fully described in Marigo et al. (1996, 1998), and Marigo & 
Girardi (2007), so that it will not be repeated here. The mixing- 
length parameter is assumed a - 1.68. 

As a matter of fact, it has long been known that the atmo- 
spheric opacity is critical in determining the position in the H-R 
diagram of a red-giant star (e.g. Keeley 1970; Scalo & Ulrich 
1975). We also recall that during the quiescent burning stages of 
both RGB and AGB phases of a low-mass star the stellar lumi- 
nosity is essentially controlled by the mass of the central core 
(and the chemical composition of the gas), being largely inde- 
pendent of the envelope mass. Adopting suitable core-mass lu- 
minosity relations available in the literature, for given value of 
the core mass and chemical composition, envelope integrations 
yield the effective temperature at the corresponding luminosity. 
We have repeated this procedure increasing the core mass - from 



0.2 M to 0.46 M for the RGB and from 0.5 M to 0.75 M for 
the AGB - and adopting different opacity tables for T < 10000 
K. 

The results for 1 .5 M and 2.0 M models with Z = 
0.02, X = 0.7 are shown in Figs. [8] and for the RGB and 
AGB tracks respectively. We have adopted low-T opacities 
from AF94, F05, LA09, and jESOPUS, and two reference so- 
lar compositions, i.e. GN93 and L03. In all cases the computa- 
tions with the opacities from jESOPUS and from the Wichita 
State University group are in close agreement, typically be- 
ing abs(log7^ F94 - logr* sopus ) ^ 0.001 dex (ranging from 

~ 5 K to ~ 20 K) and abs(log r F ° 5 - log r* sopus ) < 0.005 
dex (ranging from ~ 10 K to ~ 50 K). The deviations from 
the results with LA09 opacities are somewhat larger, 0.005 ^ 
abs(log r^ 09 - log rJ. S0PUS ) £ 0.02 dex (ranging from ~ 50 K 
to ~ 100 K). In this respect it should be recalled that in the T e ff- 
range considered here, 3.4 log(r e (j) S= 3.7, the main opacity 
contributors are the absorption by H and Thompson e~ scat- 
tering (the concentration of water vapour is still relatively low 
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even at the lowest temperatures; see Fig. [5]), so that differences 
in opacities are likely due to differences in the description of 
the H opacity, and/or in the density of free electrons, which in 
turn may be affected by differences in the partition functions of 
the ions with low-ionisation potentials. Anyhow, the temperature 
differences among the RGB and AGB tracks are in most cases 
lower than the current uncertainty affecting the semi-empirical 
r eff -scale of F-G-K-M giants (cr ~ 60 - 80 K; e.g. Ramirez & 
Melendez 2005; Houdashelt et al. 2000). 

4.2. Varying C-N-0 mixtures 

In several situations Rosseland mean opacities for non-scaled 
solar abundances should be used. One of these cases applies, 
for instance, to stellar models in which the surface abundances 
of C, N, and O are altered via mixing and/or wind processes. 
A remarkable example corresponds to the TP-AGB phase of 
low- and intermediate-mass stars, whose envelope composition 
may be enriched with primary carbon (and possibly oxygen) 
via the third dredge-up, or with newly synthesised nitrogen by 
hot-bottom burning. As a net consequence, the abundances of 
C, N, and O as well as their abundance ratios may be signifi- 
cantly changed compared to their pre-TP-AGB values (Wood & 
Lattanzio 2003). Most critical is the variation of the surface C/O 
ratio, which controls the chemistry of the gas at the low tem- 
peratures typical of the atmospheres of AGB stars (e.g. Marigo 
2002). 

Indeed, one of the aims of the present work is to provide 
a flexible computational tool to generate RM opacities for any 
value of combination of the C-N-0 abundances, hence C/O ratio. 

Figure [TOl shows clearly that big changes in kr are expected 
at low temperatures, say log(T) < 3.5, when passing from an O- 
rich to a C-rich chemical mixture. For instance, at log(T) = 3.3 
RM opacities of a gas with C/0= 1 .3 become much larger than in 
the case with C/0= 0.49 at lower densities, -8 ^ log(fl) 5 -3, 
while the trend is reversed at increasing density, log(R) > -3. 
This fact is extremely important for the consequences it brings 
about to the evolutionary properties of C stars (see e.g. Marigo 
& Girardi 2007; Cristallo et al. 2007; Marigo et al. 2008; Weiss 
& Ferguson 2009; Ventura & Marigo 2009). 

In this context we will analyse in detail the impact of chang- 
ing the C/O ratio in a gas mixture, thus simulating the effect of 
the third dredge-up in TP-AGB stars. 

4.2.1 . Molecular chemistry: the key role of the C/O ratio 

FigureQT]illustrates the abrupt change in the chemical equilibria 
when the C/O ratio passes from below to above unity, in a gas 
with log(r) = 3.3 and log(fl) = -3 (log(p) = -11.1). From 
a more careful inspection of Fig. [TT] we see that the abundance 
curves of the O-bearing molecules (top panel) and the C-bearing 
molecules (bottom panel) follow mirror trends, exhibiting two 
sudden changes of values at C/O ~ 0.93 and C/O ~ 1.0. We 
may say that these two C/O values bracket the transition region 
between the O-dominated and the C-dominated chemistry. As 
discussed by Ferrarotti & Gail (2002) the abrupt changes in the 
chemical equilibria at C/O « 0.93 and C/O ~ 1 .0 respectively 
correspond to the critical values of the carbon abundance 
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3.2 3.4 3.6 3.8 4 4.2 4.4 



log(T) 

Fig. 10. Comparison of RM opacities relative to two gas mix- 
tures with Z re f = 0.02, X = 0.7 but different C/O ratios, namely 
C/0= 1.3 and C/O=C/O = 0.49 according to GS98 solar com- 

C/O— 1 3 

position. The colour map shows the difference log(A- R ' ) - 

log(A: R '' ~ ' 49 ) throughout the standard location in the log(T) - 
log(7?) diagram of one opacity table computed with jESOPUS. 
The contour lines corresponds to differences A log(A- R ) multiple 
of ±0.25 dex. Note the large deviations occurring in the low-T 
region dominated by molecular absorption. 



The existence of e^ 1 "' 1 and £ c nt ' 2 can be understood consider- 
ing the extraordinary high bond energies of the two monox- 
ide molecules CO and SiO, i.e. £ B (CO) = 11.16 eV and 
£B(SiO) = 8.29 eV, as well as the usually large concentra- 
tions of the involved species, i.e. C, O, and to a less extent Si. 
Following Ferrarotti & Gail (2002) for temperatures T $ 1 500 
K, at which dust is expected to condensate, one must also con- 
sider the contribution of another strongly-bound molecule, SiS ( 
£B(SiS) = 6.46 eV), so that the first critical carbon abundance 
should be redefined as s^ 1 "' 1 = eo ~ £si + Since this study 
deals with the gas chemistry for log(T) > 3.2 (i.e. without dust 
formation) in the following we limit our discussion to the case 
described by Eq. (|35l l. 

In most cases the bond strength of CO mostly determines 
the chemical equilibria: as long as ec < fi™'' , the excess of oxy- 
gen atoms, so - sq, is available for the formation of O-bearing 
molecules - such as SiO, H2O, TiO, VO, YO, etc. -, while as 
soon as sc > s c nt ' 2 , i e- C/O > 1, the situation is reversed and 
the excess of of carbon atoms, ec ~ «o, takes part in C-bearing 
molecules such as CN, HCN, C2, C2H2 , SiC, etc. This also ex- 
plains why, unlike the others, the abundances of the molecules 
involving the carbon monoxide, like CO itself and HCO, show a 
flat behaviour with the C/O ratio. 

The situation is somewhat different in the transition interval, 
gp 111 ' 1 £ ec S= 4'"' 2 ' where the molecular pattern is controlled 
also by SiO, in addition to CO. The C, O, and Si atoms are 
now almost completely absorbed in the CO and SiO monoxides, 
which are the most abundant molecules, as shown in Fig.[TT] In 
other words, the excess of oxygen atoms over carbon is trapped 
in the molecular bond with silicon, which accounts for the first 
abundance drop of the other O-bearing molecules at C/O 0.93. 

It is clear from Eq. ( |35l l that the value of (C/0) cr i ti i de- 
pends on the assumed oxygen and silicon abundances. In prin- 
ciple any change in the ratio Si/O would correspond to a differ- 
ent (C/0) cr it.i. As a reference case, it is instructive to compare 
the results for different choices of the solar abundances. They 
are listed in Table [3] Passing from the AG89 to the most re- 
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Fig. 12. The same as in Fig. [4] but for gas mixtures with C/O=0.97 (upper panels) and C/O=1.30 (bottom panels) and log(T) = 3.3. 
Note that in the spectral range relevant for the RM opacity, the total monochromatic coefficient is affected by heterogeneous sources 
(e.g. TiO, ZrO, CO, H2O, CN, Thomson e~ scattering) for C/O=0.97, while absorption by C-bearing molecules dominate for 
C/O=1.30. 



cent GAS07 compilation, the C/O cl ; t; i decreases from ^ 0.96 
to ^ 0.93, implying that the transition from the O- to the C- 
dominated chemistry takes place over a wider range of the C/O 
ratio, i.e. ~ 0.93 - 1 for GAS07 in place of ~ 0.96 - 1.00 for 
AG89. As we will see later in this section, the knowledge of this 
critical ratio is of crucial importance since it defines the onset of 
the transition between two chemical regimes, with consequent 
dramatic effects on the corresponding RM opacities of the gas 
(see for instance Figs. [TBI and [TBI). 



4.2.2. Opacity sources at increasing C/O ratio 

The extreme sensitiveness of the molecular chemistry - for 
log(T) 55 3.5 - 3.6 depending on the density - to the C/O pa- 
rameter has striking consequences on the low-temperature gas 
opacities, as shown in Fig.[T2l relative to \og(T) = 3.3 and three 
values of the R parameter. This figure can be interestingly com- 
pared with the bottom panels of Fig. |4j describing the case of 
an oxygen-rich scaled-solar chemistry. For instance, we see that 
at log(T) = 3.3 and \og(R) = -3 (bottom-mid panel of Fig.fTZll 
the total monochromatic coefficient k{v) for C/O = 1 .3 is mostly 
determined by the absorption bands of molecules such as HCN 
and CN, while in a gas with the same thermodynamic conditions 
and solar C/O - 0.5, the dominating species are H2O, TiO, and 
VO (see bottom-mid panel of Fig. @J. 



At the same temperature and density, and for C/O = 0.97 
(upper-mid panel of Fig.fTZli the total coefficient k{v) is, on aver- 
age, lower than in the other two cases, being mostly affected by 
the absorption bands of CO, while the gaps in between are pop- 
ulated by the weaker molecular bands of H2O, SiO, ZrO, TiO, 
etc. At lower densities (\og{R) - -8; upper-left panel of Fig.lTZb 
Rayleigh scattering from neutral H and Thomson scattering from 
free electrons fill the spectral intervals between the CO absorp- 
tion bands, while at higher densities {\og(R) = 1; upper-right 
panel of Fig. [T2| > the total monochromatic coefficient is com- 
pletely dominated by molecular absorption, with a sizable con- 
tribution by CIA(H2/H2) at A ^ 2 /mi, just in correspondence of 
the maximum /l max of the weighting function of the Rosseland 
mean (see Eq.[28l. 

The sharp changes in the chemistry and monochromatic co- 
efficient k{v) as a function of C/O impact as much strongly on 
the integrated RM opacity kr, which is evident in Figs. |T3l— |T6| 

For the same two C/O values considered above, Fig. Qj] 
shows the contributions of different opacity sources to the RM 
opacity as a function of the temperature (and assuming \og{R) = 
-8, -3, 1). An instructive comparison with the results for a 
scaled-solar chemistry can be done with the help of Fig. In the 
case with C/O = 0.97 (upper panels of Fig. [TBI Rayleigh scat- 
tering from hydrogen and Thomson scattering from free elec- 
trons dominate for \og(R) = -8, becoming comparable with the 
molecular sources for \og(R) = -3. Moreover, we notice that 
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Fig. 13. The same as in Fig. [5] but for gas mixtures with C/O=0.97 (upper panels) and C/O =1.3 (bottom panels), and zoomed 
into the molecule-dominated temperature region. Note the various opacity bumps of the C-bearing molecules in the C/0= 1.30 case, 
while comparable contributions from both O-rich and C-rich molecules are present in the C/O=0.97 case. 



at this C/O value, representing the transition between different 
chemistry regimes, the opacity pattern is quite heterogeneous as 
it includes the contributions from both O-bearing and C-bearing 
molecules. For instance, we see that H2O is important at lower 
temperatures, CN shows up at larger temperatures, while CO 
contributes over a larger temperature interval. 

In the case with C/O =1.3 (bottom panels of Fig. [T3l the 
most noticeable features at different densities are the following. 
At \og(R) = -8 and log(r) ^3.3 the largest contribution come 
from C3 (and CN, C2), while at larger temperatures the elec- 
tron scattering dominates. At \og(R) = -3 the high and broad 
opacity bump of CN that dominates the RM opacity over a wide 
temperature interval, 3.30 ^ log(T) ^ 3.55, while the C2H2 con- 
tribution is prominent for log(T) ^ 3.30. In addition, other C- 
bearing molecules (C2, C3, HCN, CO) provide non-negligible 
contributions to the RM opacity. Finally, at \og(R) = 1 the poly- 
atomic molecule C2H2 is the most efficient contributor to kr for 
log(T) < 3.4, while the hydrogen anion becomes prominent at 
higher temperatures. 

The complex behaviour of the RM opacities as a function 
of the C/O ratio is exemplified with the aid of Fig. [14] for 
3.2 ^ log(T) ^ 3.6, the temperature range in which molecules 
become the most efficient radiation absorbers. It turns out that 
while the C/O ratio increases from 0.1 to 0.9 the opacity bump 
peaking at (log(T) ^3.3 for \og(R) = 3) - mostly due to H2O 
- becomes more and more depressed because of the smaller 
availability of O atoms. Then, passing from C/O = 0.9 down 



to C/O = 0.95 the H2O feature actually disappears and kr dras- 
tically drops by more than two orders of magnitude. In fact, at 
this C/O value the chemistry enters the transition region already 
discussed (see Fig.fTTTi. so that most of both O and C atoms are 
trapped in the CO molecule at the expense of the other molecu- 
lar species, belonging to both the O- and C-bearing groups. At 
C/O = 1 the RM opacity increases at the lowest temperatures, 
log(T) S 3.3 , while a sudden upturn is expected as soon as C/O 
slightly exceeds unity, as displayed by the curve for C/O = 1 .05 
in Fig. [14] This fact reflects the drastic change in the molecu- 
lar equilibria from the O- to the C-dominated regime. Then, at 
increasing C/O (1.1, 1.2, 1.5, and 2.0) the opacity curves move 
upward following a more gradual trend, which is related with the 
strengthening of the C-bearing molecular absorption bands. 

An enlightening picture of the dependence of the RM opacity 
on the C/O ratio is provided by Fig. [15] which displays the map 
of log(/CR) at varying temperature and C/O, for fixed log(^) = 
-3. In this diagram the drop in opacity marking the transition 
region between the O-rich and C-dominated opacity is neatly 
visible as a narrow vertical strip of width 0.95 ^ C/O 55 1.00 
(assuming GS98 as reference solar mixture) for temperatures 
3.2 < log(r) 55 3.35. This C/O range exactly coincides with 
the transition interval, (C/0) cr i t ,i ^ C/O ^ (C/0) cr i tj 2, between 
the O- and C-dominated chemistry. As already mentioned, the 
lower limit C/O cl -i t ,i is particularly sensitive to the abundance of 
silicon relative to oxygen. In respect to this, Fig. [16] shows an 
enlargement of the opacity map over a narrow interval around 
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Fig. 11. Concentrations of several gas species as a function of 
the C/O ratio, in a gas mixture with log(T) = 3.3, log(7?) = -3 
(or equivalently log(p) = -11.1), Z re f = 0.02, and X = 0.7, and 
adopting the GAS07 reference solar partitions. The increase of 
C/O follows that of C, while O abundance is kept unchanged. 
The actual metallicity Z also increases with C. The molecules 
are divided into two groups, namely: O-bearing molecules (top 
panel) and C-bearing molecules (bottom panel). Note the sharp 
change in molecular concentrations at C/O ~ 1 . 



C/O = 1, for two choices of the reference solar composition, 
i.e. AG89 and GAS07. It is evident that the opacity dip affects 
a larger C/O range in the case of GAS07 as it corresponds to 
a higher ratio, (Si/O) = 7.079 10~ 2 , compared to AG89 with 
(Si/O) = 4.168 10~ 2 . Once chosen the reference solar mixture, 
one should take this feature into account when computing RM 
opacity tables at varying C/O ratio, in order to have a good sam- 





11 

" 0.1 


1 1 1 1 1 1 


1 l 1 1 1 1 


-2 




\2.0 












^3 


l^o>^X\\\ 


\ VY 




O 


_ V\\ V 




- 


-4 








-5 


-X95 
1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 


1 1 1 1 1 



3.2 



3.3 



3.4 3.5 
log(T) 



3.6 



3.7 



Fig. 14. Rosseland mean opacity as a function of temperature, 
assuming log(7?) = -3, and at increasing C/O, from 0.1 up to 2.0. 
The reference composition is defined by (Zref = 0.02, X = 0.7) 
and assuming the metal abundances scaled-solar to the GAS07 
mixture. The abundance of carbon is made vary (hence the actual 
Z), while keeping unchanged that of oxygen. 
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Fig. 15. Rosseland mean opacity as a function of the tempera- 
ture and increasing C/O. adopting the GAS07 solar mixture, and 
assuming Z le f = 0.02, X = 0.7, and \og(R) - -3. The abundance 
of carbon is made vary accordingly to the current C/O ratio (so 
that the actual metallicity varies as well), while that of oxygen is 
kept fixed at its scaled-solar value. 



pling of the critical region, and avoid inaccurate interpolations 
between grid points belonging to different regimes. 

Going back to Fig. [15] we also notice that in the 3.4 < 
log(T) S5 3.6 the RM opacity increases with C/O. This fact is 
due to the increasing contribution from the CN molecule, which 
is one of the relevant opacity sources in this temperature interval 
(see bottom-middle panel of Fig. [5]for C/O = 0.54, and Fig. [13] 
for C/O = 0.95 and C/O = 1.3). It is worth remarking that the 
effect on the H opacity due to the increased carbon abundance 
is quite modest and only affects the opacity for log(T) > 3.6, 
when ionised carbon is expected to provide some fraction of the 
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Fig. 16. The same as in Fig.[l5j but zoomed into a narrower inter- 
val around C/O = 1 . The reference solar compositions are AG89 
(top panel) and GAS07 (bottom panel). 
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Fig. 17. Rosseland mean opacity as a function of the temperature 
in a gas with Z re f = 0.02, X = 0.7, and log/? = -3. The adopted 
chemical composition is characterised by C/O = 1.3, for vari- 
ous compilations of the reference solar mixture, namely: GN93, 
GS98, H01, L03, GAS07,and C09. In each case the actual metal- 
licity Z > Z re f because of the increase in C abundance. 
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Fig. 18. Comparison between our RM opacities and the data 
from Lederer & Aringer (2009), in terms of log(/c™ thor ) - 
log(/<^ sopus ). A few Contour lines are plotted with the corre- 
sponding values (in dex). The chemical mixture is defined by 
Zref = 0.02, X = 0.7, and C/O - 1.49. The reference solar com- 
position is L03. 
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Fig. 19. Relation between the total metallicity Z = fz Z re {, the 
C/O ratio and the degree of a-enhancement [ar/Fe], for the three 
chemical mixtures A, B, and C defined in Sect. 14.31 and adopting 
Zref = 0.02. The reference solar mixture is GS98. 



available free electrons (see Fig.l22l. A more exhaustive consid- 
eration of this point is given in Sect. 14.31 when discussing the 
case of a-enhanced mixtures. For larger temperatures the dif- 
ferences in opacity at increasing C/O progressively reduce and 
practically vanish for log(T) > 3.7, when the opacity is con- 
trolled by the hydrogen bound-free and free-free transitions. 

Let us now briefly comment the sensitiveness of the results 
to the reference solar mixture. To this aim Fig. [T7] illustrates 
the trend of RM opacity as a function of the temperature in a 
carbon-rich gas (C/O = 1.3) with the same Zef = 0.02, but dif- 
ferent choices of the solar composition. The differences show 
up for log(T) 55 3.65 and in most cases are modest, thus con- 
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firming the key role of the C/O ratio in determining the basic 
features of the molecular opacities. Another point which de- 
serves some attention is the behaviour of the RM opacity in 
the 3.55 ^ log(r) ^ 3.65 interval, which is affected mainly by 
the CN molecular bands and the negative hydrogen ion H . A 
detailed discussion of this point has been already developed in 
Sect. ED 

4.2.3. Practical hints on interpolation 

At given metallicity Z and partitions of the metal species XrfZ, 
interpolation between pre-computed opacity tables is usually 
performed as a function of the state variables (e.g. T and R) and 
the hydrogen abundance X. 

When dealing with chemical mixtures with changing ele- 
mental abundances, as in the case of the atmospheres of TP- AGB 
stars, one has to introduce additional independent parameters, in 
principle as many as the varying chemical species. 

Let us consider here the most interesting application, that is 
the case of TP- AGB stars which experience significant changes 
in the surface abundances of CNO elements, hence in the C/O 
ratio. Suppose, for simplicity, to have a chemical mixture with 
C/0> 1. Correct interpolation requires that not only the car- 
bon abundance Xq is adopted as independent parameter, but 
also the C/O ratio given its crucial role in the molecular chem- 
istry and opacity (see Figs. [TT1and[T4li. In addition, one should 
pay attention to the drastic changes in kr in the proximity of 
C/0= 1. The narrow opacity dip, delimited by the boundaries 
C/O cl it,i = 1 - esi/eo and C/O c rit,2 = 1 (see Figs. fTBTTBT l. should 
be sampled with at least 1 or 2 opacity tables, to avoid substantial 
mistakes in the interpolated values. 

A useful example of an interpolation scheme suitable to treat 
the complex chemical evolution predicted at the surface of TP- 
AGB stars undergoing both the third dredge-up and hot-bottom 
burning can be found in Ventura & Marigo (2009), where the 
grid of pre-computed opacity tables covers wide ranges of C-N- 
O abundances (and C/O ratio). Following the formalism intro- 
duced in Sect. [3] the adopted independent parameters (besides 
T, R and X) are the variation factors fc, fc/o, and (defined by 
Eq. |3T1 >, which are assigned values both > 1 (i.e. enhancement) 
and < 1 (i.e. depletion) to account for the composite effect on the 
surface composition produced by the third dredge-up and hot- 
bottom burning. In fact, the C/O ratio may initially increase due 
to the the third dredge-up and then decrease when hot-bottom 
burning consumes carbon in favour of nitrogen. 

Finally it should be remarked that, when dealing with C- 
rich mixtures, adopting both fc and fc/o (rather than either fc 
or fc/o) as independent parameters allows more robust results, 
since the interpolation is piloted by both the actual carbon abun- 
dance (mainly affecting the strength of the opacity curves) and 
the actual C/O ratio (mainly influencing the morphology of the 
opacity curves; see Fig. PHI. 

4.2.4. Comparison with other authors 

Finally, we close our discussion on the RM opacities for C-rich 
mixtures by comparing our results with the data calculated by 
Lederer & Aringer (2009). Figure[18]shows an example for a gas 
mixture characterised by Z re f = 0.02, X = 0.7, and C/O = 1.49. 
In general, the agreement between the two calculations is rea- 
sonably good, but worse than that for scaled-solar mixtures (see 
Fig-El panel b). The largest differences show up at the lower tem- 
peratures, where the RM opacity is dominated by the CN, C2H2, 



C2, HCN, C3 molecular bands. This migth appear a bit odd since 
both sets of calculations adopt basically the same molecular data 
(see Table |2). 

In the range 3.2 H log(T) ^ 3.4, compared to Lederer & 
Aringer (2009), jESOPUS predicts larger RM opacities (up to 
0.1/0.2 dex) across a strip with -7 ^ log(/?) 5 3, and lower 
values (up to =s 0.3 dex) for \og(R) > -3. 

One likely motivation of the former difference is that the 
scaling introduced by LA09 to the original gf values in the 
C2 line list (Querci et al. 1974) is not included in our calcula- 
tions. As discussed by LA09 (see their figure 10) not applying 
this correction to the line strengths of C2 causes an increase of 
log(ATR) up to a 0.1 dex, which is just what we get in terms of 
log(A^ A09 ) - log(/f^ SOPUS ) in that particular region of the dia- 
gram. On the other hand, more recently Aringer et al. (2009) 
have shown that omitting this scaling modification to the origi- 
nal C2 line list improves the comparison between synthetic and 
observed colours of carbon stars (see their figure 15). 

The latter discrepancy between LA09 and iESOPUS at larger 
densities has not a clear reason at present. We note that in this 
region of the log(T) - log(R) diagram, the dominant contribu- 
tion to the RM opacity is provided by C2H2 (see bottom pan- 
els of Fig. [T2l . We are currently investigating possible differ- 
ences among the partition function and/or dissociation energy 
of this molecule, adopted in the EOS calculations by LA09 and 
yESOPUS. 

4.3. a-enhanced mixtures 

We will analyse a few important aspects related to RM opacities 
of a-enhanced mixtures, i.e. characterised by having [a/Fe] > 0, 
according to the notation (in dex): 

[a /Fe] = log&)-log(|^) (36) 

where X a<e and Xp e are the total mass fractions of the a- 
elements and Fe-group elements, respectively. In the following 
we allocate O, Ne, Mg, Si, S, Ca, and Ti in the a-group, while 
V, Cr, Mn, Fe, Co, Ni, Cu, and Zn are assigned to the Fe-group. 
It should be noticed that, since Fe is by far the most abundant 
element of its group, the ratio [a/Fe] calculated with Eq. d36t 
coincides with the ratio computed using the abundances in num- 
ber fraction: 

[a/Fe] = log log [^-) (37) 

For simplicity in our discussion we take as selected elements all 
a-elements which are given the same [a/Fe] > 0. However, it 
should be remarked that any other prescription, concerning both 
the selected elements and the corresponding [X, /Fe] (i.e. positive 
or negative), can be set by the user via the iESOPUS interactive 
web page. 

First of all, we call attention to the fact that a given value of 
the ratio [a/Fe] is not sufficient to specify the chemical mixture 
unambiguously. The same degree of a-enhancement may corre- 
spond to quite different situations, as exemplified in the follow- 
ing. 

Adopting the formalism introduced in Sect. [3] and introduc- 
ing the quantity fz = Z/Z re f (gz = sz/sz, cr ), we define three 
different a-enhanced compositions that, in our opinion, may de- 
scribe possibly frequent applications. They are characterised as 
follows (considering the metal abundances expressed in mass 
fractions): 
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- Mixture A: Z — Zef hence fz — 1; fi > for ar-elements 
(enhanced group); fi<0 for any other element (depressed 
group). In this case the fixed group (with fi = 0) is empty. 

- Mixture B:Z — Z le f hence fz - 1; /; > for ar-elements (en- 
hanced group); fj < for the Fe-group elements (depressed 
group); fi = for any other element (fixed group). 

- Mixture C. Z > Z re f hence fz > I; fi > for ar-elements 
(enhanced group); fi<0 for any other element (depressed 
group). In this case the fixed group (with /; = 0) is empty, as 
for mixture A. 

For each of the three mixtures considered here, Table [4] lists 
the variation factors, fi and g,-, of the most relevant elements, i.e. 
C, N, O, Fe-group elements, and the metallicity parameter, fz 
and gz, as a function of a few selected [cc/Fe] values. The gen- 
eral analytical derivation of the abundance variation factors as a 
function of the selected [X, /Fe] for the three kinds of mixtures is 
detailed in Appendix |Cl Figure [19] displays the expected trends 
of C/O and fz at increasing [ff/Fe] assuming Z le f = 0.02 and the 
GS98 solar composition. 

For a given [ar/Fe] value, the three mixtures have distinc- 
tive abundance features when compared to the reference com- 
position, i.e. with Z = Z le f and scaled-solar partitions of metals. 
In particular, for their relevance to the resulting RM opacity, it 
is worth considering the changes in the CNO abundances, and 
mostly in the C/O ratio. 

- Mixture A is depleted both in the iron-group elements as well 
as in carbon and nitrogen. For instance, at [ar/Fe] = 0.4, 
the abundances of the Fe-group elements are almost halved 
and the same applies to C and N, while O is augmented by 
~ 23%. As a consequence the C/O ratio decreases signifi- 
cantly, passing from (C/O) ~ 0.49 down to (C/O) ~ 0.19. 
In general, the ratio C/O lowers considerably at increasing 
[ar/Fe]. 

- Mixture B is depleted in the iron-group elements, while C 
and N are left unchanged. At [a/Fe] = 0.4, the abundances 
of the Fe-group elements are depressed by ~ 43%, while O is 
increased by only ~ 6%. In this case the C/O ratio is just little 
affected, changing from (C/O) ~ 0.49 to (C/0) ~ 0.46. In 
general, the ratio C/O slightly decreases at increasing [a/Fe]. 

- Mixture C has the same characteristics of mixture A in terms 
of metal partitions, i.e. (X,/Z re f)A = (X//Z)c, but with a dif- 
ferent metallicity. It follows that the C case shares with A 
the same elemental ratios, so that the C/O declines signif- 
icantly at increasing [a/Fe], while the total metallicity in- 
creases. For instance, at [ar/Fe] = 0.4 mixture C corresponds 
to a metallicity Z 2Z re f (see Tableland Fig. |T9b. 

The aforementioned differences in the chemistry among the 
mixtures A, B, and C affect the resulting RM opacities, as dis- 
played by Fig.|20]for = 0.02 and Fig.[2T]for Z = 0.0001, 
both assuming X = 0.7 and \og(R) = -3. Let us first discuss 
the results of the A and B cases with Z le f = 0.02. With re- 
spect to mixture A (left panel of Fig. l20l >. we see that at in- 
creasing ar-enhancement, the opacity variations show up with 
opposite trends in two temperature intervals, namely: at inter- 
mediate temperatures, 3.50 5= log(T) $ 3.65, and at lower tem- 
peratures, 3.2 ^ log(r) ^ 3.4. Specifically, the opacity knee 
at log(r) =! 3.55 slightly smooths, while the opacity bump at 
log(T) a 3.3 becomes more prominent with increasing [ar/Fe]. 

As already discussed in Sect. 14.11 in the 3.5 log(T) £ 3.6 
interval the most effective opacity source is the negative hydro- 
gen ion (see lower middle panel of Fig. |5j, which positively 
correlates with the electron density, n e . Figure [22] shows that 
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Fig. 22. Contributions of free electrons £ = n e j/n e , normalised 
to the total electron density n e , provided by different ions as a 
function of temperature, in two gas mixtures with (Z = Z re f = 
0.02, X = 0.7, log(R) = -3) and different partitions of the 
ff-elements, namely: [ar/Fe] = 0.0 (scaled-solar abundances; 
solid lines), and [c/Fe] = +0.6 (ar-enhanced mixture of type 
A; dashed lines). The arrows indicate the increasing/decreasing 
trends when passing from [a-/Fe] = 0.0 to [ar/Fe] = +0.6. The 
highest curve (in red) displays the ratio n eM ln e ^, i.e. the electron 
density of the ar-enhanced composition relative to the scaled- 
solar case. The reference solar mixture is GS98. 



in this temperature range the principal electron donors are el- 
ements with relatively low-ionisation potentials, mainly Mg, Si, 
Fe, Al, Ca, and Na, which involve both the enhanced group and 
the depressed group. For this reason, it turns out that in the a- 
enhanced mixture of type A the decreased number of electrons 
contributed by Fe (together with C, Na, Al, Cr, Ni) is practically 
counter-balanced by the increased number of electrons removed 
from the c-atoms such as Mg, Si, and S. The net effect is just 
a very little reduction in the electron density. In the case exem- 
plified in Fig.[22]even a large c-enhancement [a/Fe] = 0.6 cor- 
responds to a reduction of n e by just ~ 6% at log(T) = 3.55. 
In turn, this small variation in n e produces a minor reduction of 
the H opacity. From a careful inspections of the results we find 
that the depression of the opacity knee at log(T) 3.55 should 
be rather ascribed to the weakening of the CN molecular ab- 
sorption bands, which reflects the depression of both carbon and 
nitrogen abundances in mixture A. In fact, at these temperatures 
and \og(R = -3) the CN contribution to the RM opacity is not 
negligible (see lower middle panel of Fig. |5). 

With respect to mixture B, we note that in the same temper- 
ature range, i.e. 3.50 ^ log(T) ;$ 3.65, the variations of the RM 
opacity at increasing [ar/Fe] are smaller than for mixture A, al- 
most negligible. In fact, in mixture B the carbon and nitrogen 
abundances are left unchanged so that the opacity contribution 
from CN is not expected to vary as well. Furthermore, the same 
arguments on the electron density, discussed for mixture A, hold 
also in this case, and the FF opacity contribution is predicted to 
change just slightly. 
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Table 4. Main characteristics of the a-enhanced mixtures described in text. The variation factors of C, N, O, and Fe-group elements, 
defined by Eqs. (|33T > for abundances either in mass fraction or in mass fraction, are indicated together with the quantities fz = Z/Z re f 
and gz - £z/£z, cr - The reference solar composition is GS98. 
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Fig. 20. Rosseland mean opacity as a function of temperature and assuming log/? = -3, for a gas with (Z re f = 0.02, X = 0.7) and 
various choices of the ratio [a/Fe], as indicated. The reference solar mixture is GS98. Results are shown for three choices of the 
chemical composition, corresponding to mixtures A, B, C described in the text. 




Fig. 21. The same as in Fig. [20] but for = 0.0001 . 



Let us now consider the temperature interval 3.2 5 log(T) ^ 
3.4, which is characterised by the opacity bump due to the 
molecular absorption bands of H2O, TiO, and ZrO. We see from 
Figure [19] (bottom panel) that the opacity peak grows at in- 
creasing [a/Fe], reflecting the decrease of the C/O ratio. In fact, 
the concomitant enhancement of oxygen and the depression of 



carbon favour the chemistry of the O-bearing molecules, thus 
strengthening the opacity contributions of H2O, TiO, and ZrO 
at those temperature. The reader should refer to Sect. 14.21 for a 
broad analysis of the dependence of the low-temperature opac- 
ity on the C/O ratio. For the same reasons, in the case of mix- 
ture B the opacity bump is practically insensitive to changes in 
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Fig. 23. Differences in RM opacities between a-enhanced mix- 
tures, with [ar/Fe] = +0.4 (top panel) and [a/Fe] = +0.8 (bot- 
tom panel), and a scaled-solar composition according to GS98. 
The a-enhanced mixtures are constructed according to the A 
scheme. In all case we assume (Z re f = 0.02, X = 0.7). Contour 
lines, with an incremental step of 0.05 dex, are over-plotted to 
guide the comparison. 



[a/Fe], since the decrease of C/O ratio is just marginal, as shown 
in Fig. [191 

Figure [23] shows the differences in terms of A log(*- R ) ex- 
pected when the chemical composition of the gas is enhanced in 
a-elements, according to mixture A. The same comments already 
spent for Fig.[20](left panel) hold here. At increasing [a/Fe] neg- 
ative deviations mostly take place in the region dominated by 
the absorption of H~, while positive variations show up at lower 
temperatures, over a well-defined region in the log(T) - \og{R) 
diagram, the boundaries of which are determined by the ther- 
modynamic conditions required to form H2O efficiently (see top 
panels of Fig. [5J-. thus becoming narrower at decreasing R. 

The case of mixture C deserves different remarks. At increas- 
ing [a/Fe] the RM opacity is predicted to be larger all over the 
temperature range 3.2 ^ log(T) S= 3.75, and the variations are 
always larger than for the other two mixtures. This fact can be 
explained simply as a metallicity effect, since the global metal 
content increases with the [a/Fe] as indicated by the fz param- 
eter (see Fig. \T% . Therefore, mixture C shares with mixture A 
the same partition of metals (i.e. the same variation factors 
see Table[4]i, but their abundances are all higher, including those 
belonging to the depressed group. The net effect is systematic 
increase of the RM opacity with [a/Fe]. 

Finally, a cautionary comment is worth being made. It should 
be noticed that the while the a-elements are the same for the 
three mixtures here considered, the differences deal with i) 
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Fig. 24. Rosseland mean opacity as a function of temperature 
and assuming log R = -3, for a gas with (Z le f = 0.001, X = 0.7). 
The reference solar mixture is GS98. Results are shown for three 
chemical mixtures, namely: i) scaled-solar abundances of met- 
als; enhanced abundances of a-elements with [a/Fe] = +0.4 
(according to mixture A; see Sect. [43); peculiar chemical pattern 
characterised by additional C-N-O-Na-Mg-Al abundance varia- 
tions superimposed to the a-enhanced mixture. See the text for 
details. 



which elements are assigned to the depressed group and to the 
fixed group, and ii) the total metallicity. The results discussed 
above show clearly that this an important point which impacts 
on the resulting RM opacities. Therefore, when using RM opac- 
ity tables one should be always aware of how the underly- 
ing a-enhanced mixture has been constructed, since his/her re- 
sults may be importantly affected. This aspect has been recently 
discussed by Dotter et al. (2007). To our knowledge available 
RM opacity tables adopt a-enhanced mixtures similar to our A 
scheme (e.g. Ferguson et al. 2005 and related website of the 
Wichita State University group). 

4.4. Other peculiar mixtures: C-N-O-Na-Mg-Al abundance 
anti-correlations 

Another relevant case is suggested by the peculiar chemical 
patterns observed in stars of Galactic globular clusters (GGC), 
being characterised by striking abundance anti-correlations be- 
tween C-N and O-Na, and Mg-Al, which are in turn superim- 
posed on a typical a-enhanced mixture (e.g. Gratton et al. 2001). 
Stellar evolution models including low-temperature RM opaci- 
ties suitable for these particular compositions have been recently 
calculated (Salaris et al. 2006; Pietrinferni et al. 2009). 

Figure [24] shows an example of RM opacities computed with 
iFiSOPUS for a gas mixture which would represent the pattern 
of extreme C-N-O-Na-Mg-Al anti-correlations, as measured in 
GGC stars (Carretta et al. 2005). The adopted abundance scheme 
is the following. We start with our reference scaled-solar mix- 
ture, characterised by Z ref = 0.001, X = 0.7 and GS98 so- 
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lar composition. Then we construct a second composition with 
[a/Fe] = +0.4 following the prescriptions for mixture A (see 
Sect. 53). The C/O decreases from (C/O) 0.49 to » 0.19, 
while the total metallicity is preserved. This fact explains the 
growth of the opacity peak due to H2O at log(T) ^ 3.4 in the 
a-enhanced mixture. The reader should go back to Sect. I4.3l for 
an extensive discussion on the differences between the two RM 
opacity curves. 

Finally we perturb the second mixture and add the C-N-O- 
Na-Mg-Al anti-correlation pattern assuming the following abun- 
dance variations in dex (see Salaris et al. 2006): log(/c) = -0.6; 
logOk) = +1-8; logCfo) = -0.8; log(/ Na ) = +0.8; log(/ Mg ) = 
-0.4; and log(/Ai) = +1.0. By doing so the total metallicity al- 
most doubles, Z = 1.97 10 3 , while the ratio [Fe/H] =s -1.5 
remains the same as in the genuine a-enhanced mixture. 

The increase in metallicity is mainly due to the augmented 
N abundance, while those of C and O both drop considerably. 
The resulting C/O ratio is now 0.31, and the total £(q + £(0) is 
decreased by 83%. This fact explains that, despite of the over- 
all increase in Z, the opacity curve of the peculiar mixture lies 
systematically lower than the others in the temperature region 
dominated by the H2O bump. 

In the temperature interval 3.4 $ log(T) £ 3.6 the differ- 
ences in RM opacity among the three curves in Fig.[24]are quite 
small and should be mainly ascribed to differences in the abun- 
dances of electron donors, which in turn affect the strength of 
the H opacity. 

4.5. Metal-free mixtures 

The last important application we discuss here deals with RM 
opacities suitable for zero-metallicity gas with a primordial 
composition. Following the standard Big Bang nucleosynthe- 
sis (SBBN), the most abundant elements to be synthesised first 
were H, He, with small quantities of D and Li, and tiny (and 
negligible) traces of Be and B. In this work we assume a pri- 
mordial mixture made up of X — 0.7521, su/sh = 4.15 x 10 -10 
(ratio of abundances by number), and Y = 1 - X— Li (hence 
Z = 0), these values being predicted by the SSBN in accordance 
with the baryon-to-photon ratio as derived by the Wilkinson 
Microwave Anisotropy Probe (WMAP; Coc et al. 2004). The 
abundances of B ad Be are reasonably neglected, since SBe/fiH 
and £b/sh < 10~ 17 according to models of primordial nucle- 
osynthesis (Thomas et al. 1993, 1994). 

Figure|25]shows the predicted chemistry of a primordial gas 
as a function of the temperature and three selected values of the 
R parameter, and correspondingly Fig. [26] illustrates the relative 
contributions of the most important opacity sources to the total 
RM opacity. It is worth noticing the following points. 

At lower densities (e.g. left panels with log(R) = -8) the 
abundance of the negative hydrogen ion H~ grows very little, 
the H2 molecule does not form efficiently even at the lowest tem- 
peratures, and the concentration of H^ is negligible (reaching a 
maximum value log(n/n tot ) 20.8 at log(T) ~ 3.35). The to- 
tal RM opacity is completely dominated by scattering processes, 
namely Thomson scattering from free electrons at higher tem- 
peratures, and scattering from hydrogen atoms at lower temper- 
atures. 

At increasing densities (e.g. going from \og(R) = -3 to 
\og(R) — 1) the abundances of most relevant species like H2 , H~, 

grow higher and higher. At intermediate densities (i.e. middle 
panel of Fig. [26) we may distinguish three different temperature 
ranges, namely: 3.2 $ log(T) 5 3.6 dominated by scattering 
from H atoms, 3.6 < log(T) ^ 3.85 characterised by the contri- 



bution of H , and 3.85 5= log(r) ^ 4.5 controlled by the contin- 
uous absorption of H (bound-free and free-free transitions). Free 
electrons are provided by H + and Li + as in the previous case. 

Finally, at the highest densities (i.e. right panel of Fig. 126) 
we notice that the RM opacity in the low-temperature region 
3.2 ^ log(T) £ 3.5 is determined by collision-induced absorp- 
tions (mainly CIA due to H2-H2 collisions); the PL opacity bump 
is prominent in the range 3.2 < log(T) ^ 4.0; and continuous and 
discrete processes due to H are dominant at higher temperatures. 

It should be remarked that Thomson scattering as well as ab- 
sorption by negative ions (i.e. PL, H 2 , He ) crucially depend on 
the amount of available free electrons. By looking at the curve 
of the electron concentration (dashed line) in Fig. [25] we see 
that, among the positive ions, three are the main electron donors 
in a primordial gas, i.e. H + , Li + , and HJ . Ionisation of hydro- 
gen atoms accounts for n e at the higher temperatures down to 
log(T) ~ 3.6-3.3 depending on the density, ionised lithium prac- 
tically provides all free electrons at lower temperatures, while 
H3 contributes free electrons only over an intermediate temper- 
ature range depending on the gas density. 

Let us first consider the case of HJ . The importance of this 
ion for the electron budget of a primordial gas has been ex- 
tensively discussed by Lenzuni et al. (1991) and Harris et al. 
(2004, hereafter also H04). In this latter paper the authors have 
pointed out that the inclusion of H^, with the most recent par- 
tition function of Neale & Tennyson (1995), may increase the 
RM opacity mostly via an indirect effect on the chemistry, i.e. 
by favouring larger concentrations of H~ and, to a less extent, 
via the direct absorption by Hi. The authors have also analysed 
possible effects on the evolution of very low-mass stars of zero- 
metallicity. In jESOPUS we have included the H^ chemistry, its 
free-free opacity, while neglecting the H^ line opacity. However, 
as shown by H04, this latter provides a small contribution (few 
%) to the RM opacity in most cases, with a peak of 15% at cer- 
tain temperatures and densities. Figure [27] (left panel) displays 
the region in the log(T) - \og{R) plane which is affected by the 
H3 via its inclusion/omission in the gas chemistry. The differ- 
ences in log(tfR) are always negative along a diagonal strip in the 
log(T) - \og(R) diagram, meaning that the neglecting HJ would 
lead to underestimate the gas opacity because we omit its contri- 
bution to n e (hence weakening the H opacity and the Thomson 
electron scattering), as well as its contribution as a true absorber 
(the free-free continuum in our computations). 

The case of Li is perhaps more interesting since the pri- 
mordial abundance of this element is predicted by the SBBN 
and accurately constrained by WMAP. An extensive analysis 
on the importance of Li for the opacity of the primordial gas 
has been carried out by Mayer & Duschl (2005), to whom the 
reader should refer for a detailed discussion. Our computations 
essentially agree with the findings of Mayer & Duschl (2005). 
From the inspection of the right panel of Fig. [22]one can see that 
even a low concentration of Li notably impacts on the resulting 
RM opacity, the effect being more pronounced at lower temper- 
atures and lower densities. For log(T) = 3.2 and \og{R) = -8 
the difference in opacity is sizable, reaching a value as high as 
Alog(/fR) =i 1.6! Figure [28] helps to get a better insight of the 
role of Li: when including it in the primordial chemistry the to- 
tal monochromatic absorption coefficient rises for A > 15 /mi 
due to the increased contribution of the Thomson electron scat- 
tering. In fact, a larger amount of free electrons is provided by 
the first ionisation of lithium, as shown in Fig. [25](left panel). 

Finally, in Figs. [29ll3T| we present a few comparisons with 
recently published RM opacity data for zero-metallicity gas, 
namely: Harris et al. (2004), Mayer & Duschl (2005), and 
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Fig. 25. Concentrations of a few atomic and molecular species as a function of the temperature in a gas with primordial composition, 
adopting Z re f = 0, X = 0.7521, and a lithium abundance of su/sh = 4.15 x 10~ 10 , and assuming \og(R) = -3. The fraction of free 
electrons, n e is depicted by a dashed black line. 
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Fig. 26. The same as in Fig.|5]but for a primordial composition with Z re f = and X = 0.7521 and assuming three of the R parameter, 
as indicated. Note the prominent bump of the CIA sources, mainly due to H2-H2 collisions, at lower temperatures in the case for 
log(/f) = l. 



Ferguson et al. (2005). In general the agreement is relatively 
good, mostly comprised within ±0.2 dex, except for the large dif- 
ferences (up to - 1 .2 - 1 .4 dex) that arise in the comparison with 
H04 and F05 at lower temperatures and densities. These discrep- 
ancies should be likely ascribed to their neglecting of Li in the 
chemical mixture, since they drastically reduce when we omit Li 
from the equation of state. We are not able to find clear reasons 
to the remaining deviations for log(T) < 3.5, temperatures at 
which Rayleigh scattering from H and H2, Thomson scattering 
from electrons, and CIA are the dominant opacity contributors 
at varying density. In general, differences in the thermodynamic 
data and input physics adopted to describe the processes listed 
in Table Q]might provide a reasonable explanation. 



5. Final remarks 
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Fig. 31. The same as in Fig. [29] in terms of log(/e^ 05 ) - 
log(A^ sopus ) between our opacity results and the tabulated val- 
ues by Mayer et al. (2005) for a metal-free mixture with Z = 0, 
X = 0.7521, and e Li /e H = 4.15 x 10" 10 . 



We have developed a new tool, jESOPUS, for comput- 
ing Rosseland mean opacities of an ideal gas in the low- 
temperature regime, 3.2 < log(T) < 4.5. The access to 
yESOPUS is made public via an interactive web-interface 
(http://stev.oapd.inaf.it/aesopus), which enables the user to spec- 
ify with large freedom the input parameters, i.e. the grid of the state variables T and R, the reference solar composition, the to 
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Fig. 27. Difference A log(/Q?) = log(KR) - log(/f^ off ) between the full RM opacity of our assumed primordial composition and the 
reduced opacity obtained either leaving out the chemistry of (left panel), or assuming a Li-free mixture (right panel). A few 
contour lines, labelled with the corresponding values (in dex), are superimposed to guide the eye. 




Fig. 28. The same as in Fig. [4] but for a primordial composition with Z le f = and X = 0.7521, and assuming \og(R) = -8 and 
log(T) = 3.3. The arrows bracket the spectral range across which the weighting function of the RM decays by a factor 1/100. The 
left panel shows the results for a lithium-free mixture, whereas the right panel illustrates the case for a primordial lithium abundance, 
as predicted by the SBBN in accordance with WMAP (Coc et al. 2004). 



tal metallicity, and the abundance enhancement/depletion of all 
chemical elements, from H to U. 

The Rosseland mean gas opacities, produced with a good ac- 
curacy (comparable to that of other opacity codes), are delivered 
in a tabular form within a reasonably short time. At present, the 
typical computation time for one table at fixed chemical com- 
position, arranged with the default T - R grid, i.e. containing 
N T x N R = 67 x 19 = 1273 opacity values, is less than 50 s 
with a 2.0 GHz processor. Such a fast performance is attained 
thanks to the optimised use of the opacity sampling method 
to describe molecular line absorption, and the adoption of pre- 
tabulated absorption cross-sections for metals (from the Opacity 
Project database). In this way the line-opacity data is suitably ar- 
ranged prior to the opacity computations, a process that, if oth- 
erwise performed on-the-fly, is in principle more accurate but at 



the cost of extremely long computing times (e.g. Ferguson et al. 
2005). 

On the other hand, several tests illustrated in the paper have 
proved that our procedure, besides being fast, is as well suitable 
to produce fairly accurate Rosseland mean opacities, to which 
the very fine spectral details are not critical as they are washed 
out, by construction, in the harmonic average of the monochro- 
matic coefficient. 

First applications of ifiSOPUS opacity tables in stellar evo- 
lutionary calculations performed with the Padova code for both 
scaled-solar (Bertelli et al. 2009), and c-enhanced mixtures 
(Bressan et al., in prep.), and with the ATON code for C-N- 
O varying mixtures along the AGB (Ventura & Marigo 2009) 
have yielded promising results. In particular, we find that the 
differences in the effective temperature of giant (RGB and AGB) 
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Fig. 29. Comparison in terms of log(*-™ 4 ) - log(Kj^ sopus ) between our opacity results and the tabulated values by Harris (2004) for a 
metal-free mixture with (Z = 0, X — 0.7). A few contour lines, labelled with the corresponding values (in dex), are superimposed to 
help the comparison. Note how much the differences become significant at lower temperatures when assuming su/sh = 4. 15 x 10 -10 
in our calculations (left panel), while they drastically reduce adopting su/en = (right panel). 




models brought about by the adoption of different opacity data 
for the same chemical composition (e.g. ifiSOPUS, Ferguson et 
al. 2005, Lederer & Aringer 2009) amount to a few tens of de- 
grees, in most cases lower than (or comparable to) the typical 
uncertainty of the semi-empirical r e ff-scale of red giants. 

We wish all interested researchers may benefit from an easy 
access to the low-temperature opacity data. Feedback and sug- 
gestions are welcome. 
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Appendix A: EOS under ICE conditions: numerical 
details 

The yESOPUS code solves the equation of state assuming in- 
stantaneous chemical equilibrium by means of the Newtwon- 
Raphson technique. We consider the N e i + 2 conservation equa- 
tions (see Sect. 12.1.2b formulated in the generic form: 



fi[n a , (a=l, 
fi[n a , (a = 1, 



•Afei), N a ,n e ] 
•A/ei), N d ,n e ] 



fN d [n a , (a =!,■■■ A/ei), N a ,n e ] 
f e [n a , (a= l,---N e i), N a ,n e ] 
/tot[« ff , (<* = !,••■ Af e i), Af a , n e ] 










0, 



which depend on the N e \ + 2 unknowns, namely: the number 
density of each neutral atom n a ipha, oc = 1, . . .N e \\ the total num- 
ber density of atoms A/ a ; and the electron density n e . Then we 
calculate the jacobian matrix J of the functions f's with respect 
to each unknown, 
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In practice, because the unknown quantities are all inherently 
non-negative functions, their logarithmic forms are adopted. 
This prevents physically unrealistic estimates from occurring 
during the iteration process. 

For each chemical element a, the abundance conservation 
equation (Eq. [T3J is conveniently written in the form: 
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fa = In 



A=l 



o, 



(A.l) 



N a £ a 



where tia is the number density of particle A which ranges over 
all species, i.e. atoms, ions, and molecules; VA,a are stoichiomet- 
ric coefficients that keep track of how many times species of type 
A contributes to the conservation equation of type a. In other 
words \>A, a represents the number of atoms of the a element con- 
tained in species A. 

The charge neutrality equation (Eq. [Tol l is expressed in the 
form: 

( N lol pz .... \ 



fe=hl 



i=l r=l 



rn A + 



+ ^—i n 



= 0. 



(A.2) 



Finally, the conservation equation of the total number density 
(Eq. W7\ is rearranged in the form: 



/tot = In 



( Kl A'raol 

^ ^ v A , a {n A + n A+ + n A -) 



a=\ A=l 



AN. 



mol 



= 0, (A.3) 



where YjA n AVA,a is extended over all molecules and quantifies 
their contribution to the conservation equation of any given ele- 
ment q-; AN 1 ™ 1 = A a - (n tot - n e ) corresponds to the excess in A a 
due to molecular formation. 

In summary, after providing a first guess to the number den- 
sities, jESOPUS sets them into the system and the jacobian ma- 
trix. In general, the guess will be inaccurate so that the functions 
fi have finite values. Denoting with F and n the entire vectors of 
the values of /, and nj, we deal with the matrix equation 



f(n + 5n) = f(n) + J(5n = 0, 



(A.4) 



which corresponds to a set of A4i + 2 linear equations for the 
first-order corrections 8n. The matrix equation is solved with the 
LU decomposition method. The corrections are then added to 



the solution vector of the number densities, n 11 



n om + c5n, and 



the process is iterated until the maximum relative change in the 
densities becomes lower than a given accuracy 5 n , i.e. typically 
10~ 5 in our computations. 



Appendix B: The frequency distribution 

Computing the RM opacity with Eq. (0 requires that the to- 
tal monochromatic absorption coefficient is evaluated at a finite 
number of frequency points. In principle the more the points, the 
more accurate the results should be. However, since we are deal- 
ing with a mean quantity one can obtain still good results using 
a relatively low number of frequency points, with the advantage 
of speeding up the computations. 

In respect to this some discussion can be found in Ferguson 
et al. (2005) who integrate over 24 000 points, and Lederer & 
Aringer (2009) who adopt 5 645 points. For the present work 
we have performed further tests to get useful indications on the 
relationship between the size of the frequency distribution and 
quality of the results, in terms of accuracy (reliability) and pre- 
cision (reproducibility) of the results. 



0.04 



0.03 - 



0.02 



0.01 




wavenumber (in 1 4 cm l ) 

Fig. B.l. The histogram of the sampling frequency distribu- 
tion with n tot = 34902 sampling points, selected following the 
scheme proposed by Helling & J0rgensen (1998). This repre- 
sents the reference distribution whence smaller frequency sam- 
ples are extracted. See the text for more explanation. 



For this purpose we proceed as follows. First, we determine a 
seed frequency distribution by adopting the scheme proposed by 
Helling & J0rgensen (1998), originally designed to optimize the 
selection of frequency points in the OS method. In few words, a 
frequency distribution produces a correct spectral sampling if it 
obeys the condition Ey(T)Av = const., i.e. expressing the con- 
stancy of the normalized energy density of the Planckian, Ey(T), 
over any arbitrary interval Ay, where v [cm -1 ] is the wave- 
number. Then, the seeked optimal distribution corresponds to the 
upper envelope the entire sample of Planckian distributions eval- 
uated at different temperatures, so that we take the maximum of 
the normalized energy density E V (T) at each v. The final distri- 
bution, shown in Fig. IB. 11 is sharply peaked at lower frequencies 
and declines exponentially at longer frequencies. 

Once the seed distribution is constructed, any other fre- 
quency grid of given size is extracted from it by using a Monte- 
Carlo technique. In our work we tested a few cases adopting 
5488, 1799, 944, 510, and 149 points. Each grid is used to com- 
pute RM opacities for two chemical compositions characterized 
by: /) X — 0.7, Z = Z re f = 0.020 and scaled-solar abundances of 
metals, and //) a carbon-rich mixture with X = 0.7, Z le f = 0.02, 
Z = 0.026 and C/O = 1.5. In the latter case carbon is made 
increase relative to its scaled-solar value, producing a net incre- 
ment of the actual metallicity. 

Then, adopting as reference opacities those obtained with the 
densest frequency grid, i.e. n tot = 5488, we evaluate the differ- 
ences, log(/C4588) - log(/c„ ), for each opacity subset computed 
with a lower frequency grid (i.e. « ; = 1799, 944, 510, and 149 
points. The results are shown in Fig. IB. 21 for mixture / and 
Fig. IB.3l for mixture //. 

We see that in most cases the differences remain small, 
within ^ ±0.05 dex, over most of the log(T) - log(R) space, and 
even with the smallest frequency set the loss in accuracy, though 
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Fig. B.2. Differences in opacities between the reference frequency grid with n tot = 5488 points, and other test cases with n tot = 
1799, 944, 510, and 149 points. The adopted chemical composition is characterized by X - 0.7, Z re f = 0.02, and scaled-solar 
abundances according to GS98. 



log Ks^g - log k 1799 log 
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Fig. B.3. The same as in Fig. IB.2I but for the adopted chemical composition which is defined by X — 0.7, Z re f = 0.02, Z = 0.026, 
andC/Q=1.5. 
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larger, is not dramatic. As expected, the biggest deviations take 
place at lower temperatures where the opacity contribution from 
molecular bands is more sensitive to the frequency sampling. 

In any case, it is worth noticing that the uncertainties brought 
about by the adopted frequency distribution are comparable, if 
not lower, with the typical differences in RM opacities computed 
with different codes (see, for instance, Figs. |7l and [T8|>. 



Appendix C: Chemical mixtures with non-solar 
K/Fe] ratios: a general scheme 

Let us first consider non-scaled-solar mixtures in which the ref- 
erence metallicity is preserved, i.e. Z = Z le f. Basing on the for- 
malism introduced in Sects. |3.2| and |4~3l we conveniently divide 
the metal species (with Z, > 3) into three groups, namely: 

- The selected elements with given y,- = [X^/Xp e ] according to 
the input specification, with abundances X? = ff Xf f ; 

f f 

- The fixed elements with abundances X\ — Xi f ; 

- The balancing elements, including all the other metals, with 
abundances X h = X b f . 

We recall that the ratios [X*/Xp e ] can be freely chosen to be 
either positive or negative. According to the adopted scheme, in 
order to preserve the metallicity the abundance variation of the 
selected elements should be compensated by the total abundance 
variation of the balancing elements. It follows that, by construc- 
tion, all balancing elements share the same variation factor f b . 

Therefore, from the condition Z = Z re f, and the definition 
of [X./Xp e ] for each of the N se i selected elements, we set up a 
system of N se i + 1 equations: 



Z re f 

Xf 



yx;~yx:~yx? 



X v 



log 



Xf 
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log 



X s 
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Xfs ) V^Fe.O/ 

which can be re-formulated with the aid of Eqs. OTb and ( T34l i: 

Vi=l,...N sel 
for the unknowns f? and f b . Let us denote with 

(CI) 



f- 

10* = j- b 



Z-0 



the partitions of metals in the solar mixture. Eventually, from 
simple analytical passages we obtain the general solution: 



Z^ + Z^ 



•/> 

k,Q 
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ft = f K' : 



(C2) 



Vi= 1, 



which only depends on the specified ratios [Xf /Xp e ] of the se- 
lected elements, and the metal partitions in the reference solar 
composition. It is useful to particularize Eq. dC.2| > for the cases 
of mixtures A and B introduced in Sect. 14.31 and finally derive 
the results for mixture C. 



Mixture A 

Since the fixed group is empty, we have 2i£, s Q 
hence: 



k,B 



= 1, 



f" = 
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1 + 2(10*-!)^ 
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Mixture B 

Since the balancing elements are those belonging to the Fe- 
group, i.e. X h k & = X Fe>0 , we get: 



Z ^ q + ^ pe ' 
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(C4) 



Fe,0 



Mixture C 

Finally, we consider the case of mixture C, in which the reference 
metallicity Z le f should not be preserved, as the actual metallicity, 
Z = fzZref, follows the total abundance variation of the selected 
elements. In this case we consider the system of equations 



fz Z re f 
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Fe 



log 
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log 
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Vi = 1, 



where we only distinguish between selected and non-selected el- 
ements. From the definitions of the abundance variation factors, 
recalling that Yij^f^~ s — Z© - Zi-^fo' anc ^ a ^ ter some manipu- 
lation, we obtain the equations 



Z 3 



Ji A i,e + Zj J® 
' j 

fi 

j^non—s 



Vi = 1, 



■Nsei 



for the unknowns ff and f" "-". We notice that formally we deal 
with exactly the same equations as those for mixture A, once the 
non-selected elements are considered in place of the balancing 
elements. Hence, the seeked solution is given by Eq. (IC.3b where 
one substitutes / with f"° n ~ s , In other words, mixture A and 
mixture C share the same non-solar metal partitions (X,/Z re f )a = 
(Xj/Z)c, but their metallicity is different by a factor fz. 



